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Chapter 1 
Introduction 


Adaptive antenna arrays have time-variable, controllable array patterns 
which enables them to suppress interference from moving sources of inter- 
ference while receiving a desired signal whose source may also be moving. 
The output of an adaptive array is the sum of the weighted (in amplitude 
and phase) outputs of a number of antennas or array elements as shown 
in Figure 1.1. The pattern of the array at any time is a function of the 
weight values at that time. The weights may be controlled according to 
some closed-loop feedback control law as in the LMS array or Applebaum 
array or they may be varied according to an open-loop weight assignment 
rule as in the SMI array [1]. In any case, an array must be given informa- 
tion that allows it to distinguish the desired signal from the interference 
signals to be suppressed. For the Applebaum and SMI arrays, the infor- 
mation consists of the desired signal direction, whereas for the LMS array, 
a reference signal that is correlated with the desired signal must be de- 
rived. Given good information about the desired signal, each of the three 
array types mentioned above is known to maximize the output signal-to- 
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interference-plus-noise ratio (SINR) under steady state conditions [1]. As 
a result, when the interference power is weak compared to the noise power 
the weights adapt to maximize the signal-to- noise ratio and the interference 
is left more or less unsuppressed. 

The problem of suppressing weak interference has been recently ad- 
dressed by Gupta and Ksienski [2] who proposed a hardware modification 
of the feedback loops which control the array weights. In this modification, 
two spatially separate antennas followed by their individual amplifiers are 
used with each feedback loop to yield decorrelated noise at the inputs of 
the feedback loop correlators. It was shown that, in theory, the modified 
feedback loops can provide the required interference suppression. Ward et. 
al. ([3], [4]) built an experimental adaptive array and verified the theoretical 
analysis. 

Later, as an alternative, Gupta [5] proposed a software modification of 
the SMI array to enhance weak interference suppression. In the modified 
SMI algorithm the sample covariance matrix is redefined to subtract out 
the effects of the thermal noise. This is done by subtracting a fraction 
F of the minimum eigenvalue of the original covariance matrix from its 
diagonal elements. This step is justified by noting that, in theory, the only 
effect of thermal noise (assuming the noise is uncorrelated, zero mean, and 
has equal power <r 2 at each array element) on the covariance matrix is an 
additive <r 2 term on each diagonal element of the true covariance matrix. 
Furthermore, if the number of elements in the array i6 greater than the 
number of signals incident on the array, the minimum eigenvalue of the 
sample covariance matrix is an estimate of the noise power er 2 . Gupta [5] 
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showed that, in theory, the modified SMI algorithm provides the required 
interference suppression assuming the true covariance matrix is known. 

In normal applications, the true covariance matrix is not known a priori 
due to incomplete knowledge of the receive antenna patterns and received 
signal characteristics. As a result, samples of the received signals (i.e. snap- 
shots) are used to estimate the covariance matrix which, in turn, is used in 
the weight calculation. Dilsavor and Moses [6] have addressed the question 
of how many snapshots are needed to achieve good estimates of the mod- 
ified SMI weights. Using Monte Carlo simulations and statistical theory, 
they found that the number of snapshots required for good estimates of the 
modified SMI weights increases with the amount of diagonal subtraction 
(i.e. with the fraction F). Thus, a tradeoff was found to exist between the 
improved suppression of weak interference provided by the modified SMI 
algorithm and the number of snapshots required to achieve that improve- 
ment. 

The modified SMI algorithm has been implemented on the aforemen- 
tioned experimental array built by Ward et.al. The experimental SMI array 
is fully adaptive with either two or three elements. It can simulate narrow- 
band signal scenarios consisting of one desired signal arriving from broad- 
side and up to two interference signals arriving from arbitrary directions. 
Since the original system built by Ward was a modified Applebaum array 
in the sidelobe canceller configuration, some hardware as well as software 
alterations were necessary in implementing the modified SMI algorithm. 
The purpose of this report is to describe the implementation of the modi- 
fied SMI algorithm on the original experimental system and then to study 
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the performance of the experimental SMI system in comparison with theory. 

The rest of the report is organized as follows. Chapter 2 offers a brief 
description of the original experimental system built by Ward. Chapter 3 
describes the implementation of SMI on the original system and includes 
a description of hardware modifications, sampling procedures, the imple- 
mented algorithm, weight normalization, performance evaluation, and sys- 
tem calibration. Experiments conducted to evaluate the performance of the 
experimental system are described in Chapter 4. Chapter 5 contains our 
conclusions. Finally, the system software is included as an appendix. 
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Chapter 2 

The Original Experimental 
System 


The experimental system built by Ward, et. al. ([3], [4]) is a sidelobe 
canceller with an unweighted main element and two weighted auxiliary 
elements which employs a discrete version of the Applebaum algorithm for 
weight control. The signal scenarios which can be tested consist of a desired 
signal arriving from broadside and up to two interference signals arriving 
from arbitrary directions. 

The system is designed to simulate the array whose auxiliary element 
weights are controlled by the modified feedback loops of Figure 2.1. In 
these modified feedback loops the effect of noise on the array weights is 
decreased by reducing the correlation between the noise components of the 
two inputs to the correlator. The reduction in correlation is achieved by 
using two spatially separate antennas, each followed by its own amplifier, 
with each auxiliary element as shown in Figure 2.1. In this configuration, 
the “3-element” array actually uses 5 antennas; one for the unweighted main 
element and two for each of the weighted auxiliary elements. For purposes 
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of parameter control and system performance evaluation, the experimental 
system can simulate the signals that would be received by the five antennas. 
Desired, interference, and noise signals are bench generated and combined 
in an array simulator to simulate the signals received by the five antennas. 
The experimental adaptive processor can also be used with signals received 
from geosynchronous satellites. In this case, the signals received by actual 
antenna elements are downconverted to 69MHz, the center frequency of the 
adaptive processor, and then fed to the experimental system. In the results 
presented in this document, the simulated signals are used. 

A block diagram of the experimental system is given in Figure 2.2. The 
signal simulator generates a desired signal and two interference signals. In 
the array simulator, these signals are combined with each other and with 
noise to form the signals that would be received at the 5 antennas. The 
array processor applies weights to the auxiliary element signals and forms 
(by summing) the adapted array output. The system operates at 69 MHz 
with a bandwidth of 6 MHz. A digital computer (PDP 11/23) is used to im- 
plement the weight control algorithm, control the various components, and 
evaluate the system performance. A description of the individual system 
blocks is given below. 

2.1 Array Simulator 

Figure 2.3 shows a detailed block diagram of the array simulator. In the 
array simulator, the incident signals are combined and thermal noise is 
added to form the signals received at each array element, such that each 
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Figure 2.1: Modified feedback loop that uses separate antennas and ampli- 
fiers to decorrelate noise at the correlator inputs. 


ARRAY 

OUTPUT 
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Figure 2.2: Block diagram of the experimental system. 
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signal contains a component due to the desired signal, components due 
to one or both interfering signals, and additive thermal noise. Thus the 
array simulator has three inputs for the three incident signals, and five 
outputs corresponding to the five antennas of the array. These five outputs 
are designated MAIN, AUX-1 SIGNAL, AUX-1 CORRELATOR, AUX-2 
SIGNAL, AND AUX-2 CORRELATOR. The MAIN output is the signal 
received by the main element of our array. 

The other outputs are the signals received by the two auxiliary elements 
of our sidelobe canceller with modified feedback loops. The designations 
SIGNAL and CORRELATOR specify the two branches of the modified 
feedback loop (Figure 2.1) of a particular auxiliary element. The SIGNAL 
branch is weighted and proceeds to the array output summer, whereas the 
CORRELATOR branch forms the input to the correlator. Note in Figure 
2.3 that the noise components injected into the auxiliary SIGNAL branches 
and the MAIN channel are all from different noise sources, and thus are 
uncorrelated. Furthermore, the noise components in the auxiliary COR- 
RELATOR branches originate from another noise source, and are therefore 
uncorrelated with the noise components of the SIGNAL branches. In Fig- 
ure 2.3 the A s are zero-phase power dividers connected as summers. The 
a’s are variable attenuators and the <£’s denote variable phase shifters. Nj 
through N4 are the noise sources. The phase shifters simulate variations of 
the interfering signal directions of arrival by varying the interelement phase 
shifts between interfering signal components of different array elements . 
There are no phase shifters associated with the desired signal because it is 
assumed to arrive from broadside and thus is received with the same phase 
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Figure 2.3: Array simulator; detailed block diagram. 









at each array element. Variable attenuators are used to control the amount 
of each incident signal received at each output channel. This is analogous 
to varying the gains of the main and auxiliary elements in the directions of 
the incident signals. Once the desired signal scenario is set, the array simu- 
lator outputs are fed to the array processor, where the signals are sampled 
and the auxiliary element weights are determined by digital computer and 
are implemented in analog. The array processor is discussed next. 

2.2 Array Processor 

Together with the system computer, the array processor forms the two 
modified feedback loops of the sidelobe canceller, through which the weight 
control algorithms are implemented. A detailed block diagram of the array 
processor is shown in Figure 2.4. Note that the auxiliary channel correlator 
branch signals are down converted to baseband and quadrature detected by 
the vector demodulators (VDMs), as is the array output. These baseband 
voltages are simultaneously sampled, analog-to-digital (A/D) converted and 
read by the system computer which implements the weight control equa- 
tion and calculates the array weights. The new weights are then (D/A) 
converted and applied to the auxiliary element SIGNAL branches as in- 
phase (I) and quadrature (Q) control voltages by the two vector modula- 
tors (VMODs). The weighted auxiliary elements are then summed with 
the main channel signal to form the array output. 

In the array processor, the I and Q outputs of each vector demodulator 
are processed prior to being sampled. A low pass filter first removes the 
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Figure 2.4: Array Processor; detailed block diagram. 




















second harmonic. The resultant baseband signals are then amplified to 
utilize the full dynamic range of the A/D converter. Track and hold devices 
allow the multiplexing of all six VDM outputs to the single A/D converter 
, so that they can be sampled simultaneously. This process preserves the 
signal correlation and noise correlation (if any) between the samples of 
different channels. For the experiments to be described here, the inputs to 
the processor are the signals from the array simulator. However, this same 
array processor may be used with signals from an actual antenna array, 
after downconversion to the array processor center frequency of 69 MHz 
[7]. The signals provided to the array simulator are described next. 

2.3 Signal Simulator 

The signal simulator synthesizes the desired signal and the two interfer- 
ing signals, which are then combined in the array simulator to form the 
received signals at each array element. In order to measure adaptive ar- 
ray performance characteristics such as interference suppression, output 
interference-to-noise ratio (INR), and output signal-to-interference-plus- 
noise ratio (SINR), it is necessary to measure separately the desired signal 
power, the interference signal power, and the noise power present in each 
of the array elements and in the array output. Pulse modulated sinusoids 
are used as the desired signal and the interfering signals to accomplish this 
objective. As shown in Figure 2.5, the modulation on one interfering signal 
is staggered from the modulation on the other interfering signal, and from 
the desired signal modulation, such that the signal occupies a different por- 
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tion of the pulse repetition period. There is also a portion of the pulse 
period when only noise (no signal) is present. The desired and interfering 
signals are therefore all uncorrelated with each other (for all interelement 
time delays of interest). The track/hold devices of the array processor are 
triggered in synchronism with this modulation envelope. Because of A/D 
conversion speed limitations, successive samples are not taken in the same 
waveform period. Instead, each sample is from a different pulse repetition 
period , but is separated by only a small time interval from the point on 
the waveform at which the previous sample was taken. Thus, an effective 
sampling rate much higher than the sampling rate possible in real time is 
achieved. By varying the delay from the start of a period to a sampling 
instant, a sequence of samples covering the entire waveform is provided to 
the system computer. The levels of a particular signal component (desired, 
interference, or noise) are obtained by signal averaging over that portion 
of the pulse repetition period corresponding to the signal component under 
measurement. In this manner, the sampled data and the computer are also 
used for steady state adaptive array performance evaluation. However, as 
discussed in the next chapter, this pulse modulation scheme is exploited 
solely for performance evaluation and not in determining the auxiliary ele- 
ment weights. The experimental modified SMI system is described next. 
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Figure 2.5: Staggered Pulse modulation of the desired and two interference 
signals. 
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Chapter 3 

Experimental Modified SMI 
System 

A number of issues were addressed in implementing the two- or three- 
element fully adaptive modified SMI array on the original system. These 
issues as well as the modified SMI algorithm itself are described in this 

chapter. 

3.1 A Hardware Change in the Original Sys- 
tem 

The modified SMI system represents a software solution to the problem 
of weak interference suppression, whereas the original system embodies a 
hardware solution (modified feedback loops) to the same problem. Orig- 
inally two antennas and amplifiers were used with each feedback loop in 
order to decorrelate the noise components of the two signals at each feed- 
back loop correlator input. The original system modeled this situation by 
using different noise sources to generate the noise components. Separate 
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noise sources for the SIGNAL and CORRELATOR branches of each ele- 
ment of the array processor are not needed in the SMI system since this 
system does not need two separate antennas for each auxiliary element. 
Instead, the noise correlation will be reduced in the software. As a re- 
sult, only three of the five outputs of the array simulator are needed by 
the three-element SMI array processor. One of these three is the MAIN 
element output of the array simulator which is fed directly to the MAIN 
SIGNAL INPUT of the array processor. The other two signals required by 
the array processor may be chosen from the remaining four array simula- 
tor outputs. Each of the two chosen signals is split into two parts using 
a two-way power divider. One output of the power divider is fed to the 
correlator branch while the other is fed to the signal branch of an auxiliary 
branch. For example, if the AUX-1 SIGNAL output of the array simulator 
is chosen to be the signal received by the second auxiliary element of our 
array then the AUX-1 SIGNAL output is split into two parts with a two- 
way power divider. One output of the power divider is fed to the AUX-2 
SIGNAL input of the array processor and the other output of the power 
divider is fed to the AUX-2 CORR input of the array processor. Care has 
been taken not to choose the AUX-1 CORR and AUX-2 CORR outputs of 
the array simulator to be received by the array processor auxiliary elements 
since to do so would create a scenario in which the auxiliary element signals 
have correlated noise components; contrary to an assumption made in the 
modified SMI algorithm derivation. To see this, recall that a single noise 
source (N2 of Figure 2.3) is used to generate the noise components of both 
the AUX-1 CORR and AUX-2 CORR outputs of the array simulator. 
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3.2 Sampling and Weighting the Main Chan- 
nel 

The fully adaptive SMI array requires that the main channel as well as the 
auxiliary channels be sampled and weighted. Since the original system was 
a sidelobe canceller, the main channel was not equipped with a VDM for 
sampling; nor does it have a VMOD for weighting the main element signal 
(see Figure 2.4). The problem of sampling the main channel is solved 
by noticing that, if the auxiliary element weights (VMOD-1 and VMOD- 
2) are set to zero, the array output signal is equal to the main channel 
signal except for some attenuation and phase shift due to the real system 
components. Thus, the main channel is sampled by setting VMOD-1 and 
VMOD-2 to zero and sampling the output of VDM-3, the array output 
vector demodulator. The attenuation and phase shift will be offset by 
sample scaling and system calibration procedures that will be discussed 
in later sections of this chapter. In order to effectively weight the main 
channel, the weights computed by the SMI algorithm are always normalized 
so that the main element has unity weight. As a result, the main channel 
is weighted without a vector modulator. 

3.3 Scaling the Samples 

A significant difference between the SMI algorithm and the original Ap- 
plebaum algorithm is that the SMI algorithm is open-loop whereas the 
Applebaum algorithm is closed-loop. In other words, the SMI algorithm 
derives weights using only the samples of the array element signals but the 
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Applebaum algorithm uses the array output to adjust the array weights. 
The input to the SMI algorithm consists of a set of K snapshots. Each 
snapshot is a vector (of length 3 for a 3-element array or length 2 for a 
2-element array) of samples resulting from a simultaneous sampling of the 
array elements. Since the SMI algorithm derivation assumes an ideal array, 
the snapshots presented to the algorithm must be identified with samples of 
the array element signals of an ideal array. To make this identification, the 
samples taken from the A/D converter must be scaled in order to compen- 
sate for losses and phase shifts through the nonideal hardware components 
of our real array. Figure 3.1 compares the block diagrams of the ideal array 
and our experimental array. The 9dB attenuators in the auxiliary channels 
represent the inherent loss through the vector modulators when the weights 
are set equal to unity. The first 6dB attenuator at the output represents 
loss through the summer El and the second represents loss through the 
4-way power divider A6 of Figure 2.4. The 9dB gain accounts for the fact 
that A5 and A6 are set 9dB higher than A1-A4 in order to use the full 
dynamic range of the A/D converter. The absolute settings of A1-A6 are 
not important since we are concerned only with relative gains. The large 
box of Figure 3.1 encloses an ideal array which “exists within” our real 
system. The dots denote the locations where samples from the ideal array 
exist; and such samples shall be called ideal samples. The samples which 
are available at the VDM outputs shall be called the real system samples. 
As seen in the figure the auxiliary element ideal samples are simply the 
corresponding real system samples attenuated by 9dB. The main element 
ideal samples are the real system output samples taken with the weights 
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set to zero and boosted by 3dB (=6+6-9). Also note that the ideal adapted 
array output samples which will be needed later to calculate output powers 
of the adapted array are the real system output samples taken with the 
weights set to their adapted values and boosted by 3dB (=6+6-9). This 
scaling of samples is done in software as soon as the samples are retrieved 

from the A/D converter. 

The scaling of samples described above represents a coarse amplitude 
adjustment in that, for example, the inherent loss through the first auxiliary 
element weight is not exactly 9dB. Furthermore, phase shifts through the 
real devices have not yet been offset. The further fine tuning which is 
required for the SMI algorithm to perform up to its abilities takes place 
in the weight calibration procedure which is discussed at the end of the 
chapter. 

3.4 The Modified SMI Algorithm 

This section describes the steps which are taken upon retrieving the scaled 
samples from the A/D converter that results in a set of modified SMI 
weights to be implemented on the array. The strategy here will be to 
state the equation used to generate the weights, define the terms in this 
equation, and thoroughly discuss how these terms are computed. 

In the Modified SMI algorithm a vector of complex weight estimates 
W K based on K snapshots are computed using 

W k =^kS (3-1) 

where the notation • will be used throughout this discussion to indicate esti- 
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Figure 3.1: Comparison of the ideal array and our experimental system. 
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mate. For example, the algorithm yields an estimate based on K snapshots 

of an optimal weight vector. For our 3-element array, the weight vector is 

defined to be A , 

w K1 ™kii ~ JWkiq 

W K = W K2 — WK2I - j'U’K2Q (3-2) 

W K3 J [ W K3 [ - jw K3 Q 

where Wki is the estimated complex weight for the first auxiliary element, 
wj (2 is the estimated complex weight for the second auxiliary element, and 
w K3 is the estimated complex main channel weight. In order to apply a 
complex weight in our real array system each VMOD i6 provided the in- 
phase and quadrature components of the proper complex weight. Equation 
(3.2) defines the in-phase and quadrature components (subscripts / and 
q, respectively) of the complex weights. The minus sign is included in 
the definition since a VMOD input signal is split (within the VMOD) into 
two components such that the quadrature-weighted component lags the 
in-phase-weighted component by 90 degrees. 

In Equation (3.1), // is an arbitrary complex scalar which is chosen here 
so that w K3 = 1 + j0. This normalization is necessary for our array system 
since the main channel is not equipped with its own VMOD. IV is the 
(3 x 3) modified sample (estimated) covariance matrix given by 

IV = - F\ min ($K)I ( 3 * 3 ) 

and S is the steering vector. The method for generating a steering vector 
is presented at the end of this chapter. 

In Equation (3.3), the standard sample covariance matrix 4>* based on 
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K complex signal snapshots {X*}^ is the block average 

♦h4EW". (3-4) 

™ fc=l 

where H indicates Hermitian transpose. Also in Equation (3.3), F is the 
fraction between 0 and 1 which is chosen to achieve a desired level of inter- 
ference suppression, A min (4»/f) is the minimum of the 3 real positive eigen- 
values of the positive definite Hermitian matrix and I is the (3 x 3) 
identity matrix. 

The procedure for computing the array weights given K complex (3x1) 
signal snapshots is clear. First, the standard sample covariance matrix 4* 
is computed as in Equation (3.4) and its minimum eigenvalue is found. 
These quantities are used in Equation (3.3) along with the chosen value of 
F to yield the modified sample covariance matrix fjf. Next, f* and the 
approximate steering vector 5 are substituted into Equation (3.1) to yield 
the complex weights Wk . Remember that fx is chosen to make wk 3 = 1 +j0. 
Finally, the I and Q components of each complex weight are found from 
Equation (3.2) and are provided to the VMODs in analog. The next task 
is to describe how the complex signal snapshots X*. are generated. 

Let Xk\, x* 2 , Zfe3, be the scaled (see Section 3.3) complex samples of the 
signals received at the first auxiliary, second auxiliary, and main channel at 
time index k. Recall, for example, that 1*3 is a scaled version of a sample 
taken through VDM-3 while the auxiliary element weights are set to zero. 
The k t/l signal snapshot X*. is a vector of these samples given by 

*fci *fcl/ + jXklQ 

Xk = = Xk2l + JXk2Q • (3-5) 

Xfc3 J [ Xfc3 / + jXk3Q . 
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In the communication signal environment, the k th sample of the n 
element will contain desired, interference, and noise components, and thus 
can be written 

M 

Xkn = XknD + X X fcn/m + X knr , (3.6) 

m — 1 

for M interference signals. Using (3.6) in (3.5), the k snapshot vector is 
written in terms of its signal component vectors 

M 

X k = X kD + X X ^im + X kv . (3-7) 

m = l 

The SMI algorithm uses K such snapshots to form the covariance matrix 
estimate 4 >k which can be written using (3.7) in (3.4) as 
K Af M 

= — X(^*£ + X Xkim + x kr ,)(x kD + X Xkim + Xkr >) H 

K fc— 1 m = l m = 1 

. K M 

= T7 T.(x kB x^ B + Z x u ~x» m + 

** fc=i m=1 

, K M M 

+ i E(x„n Z x?, m + Z *«■»*£) 

^ fc— 1 m— 1 m=l 

+ -77 Yl( XkDX kv + Xk n X ko) 

** fc=l 

. K M M 

+ -uZ(X k DZ X klm X x kIm x” D ) 

ft fc-1 m= 1 m=l 

1 K M M 

+ ^Z £ ZX U ~,X« m . (3-8) 

^ fc — 1 f»| =1 m2 =1 

mi ^mj 

The first line of the final expression in (3.8) approaches the true covariance 
matrix (which we are estimating) whereas the last four lines approach zero 
as K —* oo. 
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When signals received from geosynchronous satellites are used as input 
to the array processor, a snapshot is formed as a vector of signal sam- 
ples generated by simultaneously sampling the vector demodulator VDM 
outputs of the array processor (see Figure 2.4) and applying the calibra- 
tion factors (see Sections 3.3 and 3.6). However, when the bench-generated 
signals (i.e. the outputs of the array simulator) are used as inputs to the ar- 
ray processor, care must be taken in forming a set of K snapshots 
so as not to take advantage of the special pulse-modulated signal scheme 
employed by our experimental system. Recall that the desired and inter- 
ference signals are separated in time (see Figure 2.5) so that their powers 
can be calculated and our system performance evaluated. Thus, if (when 
using bench-generated signals) a snapshot is formed by simply combining 
the samples taken at a single instant of time then no single snapshot would 
contain both desired and interference signal components. As a result, the 
last two lines of (3.8) would not be present in the resulting covariance ma- 
trix estimate. Although this lack of desired-interference crossterms would 
not change the asymptotic value of the covariance estimate it would give an 
unrealistic (overly-optimistic) picture of the performance of the array based 
on a finite number of samples. To remedy this situation, a single snapshot 
Xt containing all signal components is “built” by simply summing three 
phase-shifted single-component snapshots; making sure that one of the ad- 
dends comes from the desired signal portion of the pulse repetition period 
(call this snapshot Xp), another (X/i) from the first interference signal por- 
tion, and the third ( A/ 2 ) from the second interference portion (see Figure 
2.5). An addend from the noise-only portion of the pulse repetition period 
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is not necessary since noise is already present in the other three addends. 
Thus, a typical bench-generated snapshot is formed as 

X T = e’*'X D + e A I n + e’*'X n . (3.9) 

The phase angles <pi, <f> 2 , and <f> 3 are i.i.d. uniform random variables on 
[0, 27 tJ . The random phases are needed to decorrelate the bench-generated 
desired and interference signals. This apparent correlation between received 
signals was a byproduct of using a single 69MHz oscillator to generate 
the carrier for the desired and interference pulse signals. The summing 
procedure (3.9) for forming bench-generated snapshots was used in [6] in 
order to study the performance of the system as a function of the number 
of snapshots. 

This report is specifically concerned with the steady state performance 
of the modified SMI system. Thus, it is to our advantage here to omit 
the last two lines of (3.8) from the covariance estimate by not using the 
snapshot summing procedure of (3.9). In other words, in this report, a 
snapshot is formed as a vector of (calibrated) samples which were taken at 
a single instant of time. In doing so we shorten the sample-taking process. 

The above presentation of the SMI algorithm was formulated in terms 
of a three-element array. The experimental array is also capable of imple- 
menting a two-element array by simply ignoring one of the two auxiliary 
elements. For this two-element configuration the algorithm must be ad- 
justed accordingly. In particular, the estimated weight vector Wk becomes 
length 2, with the first entry being the auxiliary element weight and the 
second entry being the main element weight. As before, the weight vector i6 
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normalized so that the main element weight (now Wkz) is unity. Similarly, 
the snapshot vector becomes length 2 with the auxiliary element sample oc- 
cupying the first entry and the main element sample occupying the second 
entry. 

It may have occurred to the reader that the calculation of input and 
output signal and noise powers must be thought out carefully so that the 
performance results obtained here relate directly to results found in the 
literature. The power calculation procedure is discussed next, followed by 
the topics of system calibration and steering vector generation to conclude 
the chapter. 

3.5 Power Calculations 

In order to quantitatively examine the steady-state performance of the ex- 
perimental modified SMI array it is necessary to calculate (from a set of 
samples) the desired and interference signal powers and the noise power in 
the adapted array output signal and the two or three input signals (depend- 
ing on whether a two- or three-element array is being implemented). Using 
these power measurements, quantities such as interference suppression (IS) 
and output signal-to-noise ratio (SNR) can be found and compared with 
those of theory. Comparison of theoretical and experimental system results 
for several steady-state experiments appear in the next chapter. This sec- 
tion explains the procedure for calculating the various signal powers, IS, 
SNR, etc. 

Figure 3.2 shows the I and Q outputs of one vector demodulator (VDM). 
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The D pulses are the desired signal, and II and 12 are the interfering signals. 
Due to the pulse modulation, the desired signal amplitudes are calculated 
from the sampled data over time interval r l5 and the interfering signal 
amplitudes are calculated from the samples over r 2 and r 3 . Samples over 
interval r 4 are used to compute the noise power in each channel and, in 
addition, are used to calculate and correct for any DC offset voltages due 
to the detector amplifiers (A1-A6 of Figure 2.4). Even though these offset 
voltages have been nulled by adjustments at the amplifiers themselves, this 
correction is done in software as a precaution against any drift that may 
occur during the course of an experiment. Left uncorrected, such offset 
voltages would cause errors in the estimated covariance matrix and thus 
degrade adaptive array performance. 

The signal level and noise power calculations are based on finite numbers 
of samples of a noisy signal and thus are estimates of the actual levels. In 
the experiments to be conducted, the suppressed interference components 
in the adapted array output signal will be as much as 30 dB below the noise 
level, making their levels very difficult to estimate accurately. The number 
of samples required for accurate estimates of the output interference signal 
levels (as well as all other signal levels) is reduced by simply turning off 
the noise generators while samples for signal power calculation are being 
taken. Of course, the noise generators are “on” when samples for weight 

and noise power estimation are taken.. 

Let x kJ i and x kJ Q be the k th scaled I and Q samples of element J for 
J = 1, 2, and 3, as defined in Section 3.4. These samples are not the very 
same samples used to form the snapshots needed for weight computation 
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Figure 3.2: Typical I and Q VDM outputs. 

but instead are different samples of the same signals. Furthermore, newly 
define and x^q to be the k lh scaled I and Q samples of the adapted 
array output signal. Suppose that this new data has been taken for k 
ranging from 1 to fH in each of the 4 time intervals r; of Figure 3.2. Then 
compute 

1 K1 

Ojp = ■— Y.XUP on r 4 (3.10) 

A 1 k=i 

on Ti (3.11) 


1 K1 

Djp = — Y x kJP - Ojp 

* 1 fc=i 
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Iljp = 

1 Kl 

177 XkJP ~ ® JP on T2 

** 1 


1 K1 

I2jp = 

ITT XkJP ~ ® JP ° n T3 

^ 1 k=l 

II 

1 K1 

— J2( x ^ p - 0jp ) 2 on T4 

K 1 A:=l 


(3.12) 

(3.13) 

(3.14) 


for J <E {1,2, 3 , 4 } and P € {I,Q}- Oji and Ojq are the estimates of the 
offsets in the I and Q components of array signal J. Dji and Djq are 
the estimates of the I and Q component amplitudes of the desired signal 
component in array signal J. flji, /l jQi J2 jj, an d I^JQ are analogous 
parameters for the first and second interference signals, respectively. Recall 
that the samples for the estimates of (3.10)-(3.13) are taken with the noise 
generators turned off. and a 2 JQ are the noise power estimates in the I 
and Q channels of array signal J . 

From these equations we see that the signal amplitude calculations for 
the desired and interference signals are obtained by sample averaging over 
the particular portion of the pulse repetition period corresponding to the 
signal under measurement. The variances of these signal level estimates 
will depend on the average noise power <x 2 present in the particular array 
signal being averaged and on the number of samples K 1 in the average. 
Specifically, the averages (3.10)-(3.13) are efficient maximum likelihood es- 
timates of the I and Q offsets and signal amplitudes. The variance of the 
offset estimate of (3.10) is <r 2 /Kl whereas the variance of the signal level 
estimates of (3.11 )-(3.13) is 2<r 2 /Kl. The noise power estimate of (3.14) is 
asymptotically unbiased [8]. 
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From the I and Q signal amplitude and noise power estimates in array 
signal J, the average power of the desired, interference, and noise compo- 
nents of complex array signal J are computed as follows: 


Pdj 

= d 2 ji + d% 

(3.15) 

Piu 

= Iljr + I\ JQ 

(3.16) 

Pl2J 

~ ji + P^jq 

(317) 


= 3(dj 7 + (Tjq). 

(3.18) 


The factor of 3 in (3.18) is included to account for the fact that each 
signal snapshot used in the weight estimate contains a noise component 
whereas the desired signal is present in just 1/3 of these snapshots as are 
the first and second interference signals. Because the I and Q components 
of narrowband noise are uncorrelated [9], the total average noise power in 
a particular channel is the sum of the noise powers at the I and Q outputs 
of that VDM. 

Once the power of each signal component (desired, interference, noise) 
in each array signal (first auxiliary, second auxiliary, main channel, adapted 
array output) has been computed any interesting performance measures can 
be easily found. For example, the output signal-to-noise ratio (SNR), first 
interference suppression (ISl), second interference suppression (IS2), and 
total interference suppression (1ST) are just 

P D4 

Pr) 4 


SNR 


(3.19) 



IS2 = #21 (3-21) 

* 1 23 

r CT = p n4 + ^ 24 ( 3.22) 

P/13 + Pm 

In Chapter 4, these performance measures will be used to compare the per- 
formance of the experimental system with that of an ideal modified SMI 
array operating in the same signal environment. Before the experimen- 
tal array can perform up to its capabilities, however, the system must be 
calibrated. System calibration is discussed in the next section. 


3.6 System Calibration 

The experimental system uses a software calibration procedure to correct 
for non-ideal characteristics (i.e. losses and phase shifts) of the real system 
components. Recall that the sample scaling of Section 3.3 represented a 
coarse adjustment of sample amplitudes. Phase corrections and further 
amplitude corrections to the samples are necessary to ensure proper array 
performance. 

The calibration procedure used here adapts the weights for the sim- 
ple scenario of a 2-element array receiving only a single interference signal 
(no noise, no desired signal). Since the scaled samples are not ideal, the 
estimated weights will not be ideal. However, ideally we know that the 
array should suppress the interference at the array output beneath mea- 
sureable levels. By scaling the amplitude and shifting the phase of the 
adapted auxiliary element weight, an ideal weight which completely sup- 
presses the interference is found. The amplitude scale factor A and phase 
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shift (f> required to transform the adapted auxiliary element weight to the 
ideal weight can be represented by the complex number c = Ae**. By apply- 
ing the “fine-tuning” sample adjustment factor l/c* to the scaled auxiliary 
element samples in subsequent experiments the resulting adapted auxiliary 
element weight is made equal to the ideal weight. To see this, suppose the 
fine-tuning adjustment factor for the first auxiliary element samples were 
known; call it l/c*. Then the k th fine-tuned snapshot X k FT may be written 
in terms of the k th coarse-adjusted snapshot X k cA (see Section 3.3) as 

[ 1 o' 

XkFT = q j X kCA (3.23) 

= CXkCA . (3.24) 

As a reminder, we are considering the two-element array in which the first 
auxiliary element sample occupies the first element of the snapshots and the 
main element sample occupies the second element of the snapshots. Thus, 
(3.23) implies that we are calibrating the first auxiliary element samples 
with respect to the main element samples. The estimated covariance ma- 
trix based on K of the fine-tuned snapshots may be written in terms 
of the corresponding covariance matrix 4» Cj 4 based on the coarse-adjusted 
snapshots. Using (3.24) in (3.4), we have 

$FT = ~77 X! X kFTXkFT 

* fc=l 

= J7E[CXkCA}[CX k CA) H 

K k=l 

= W^ cx 

K fc=l 
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= C 1? t X„caX“ ca 

k=l 

= C$ C aC h . 


t H 


(3.25) 


Using (3.25) in (3.3) with F = 0 (standard SMI) and substituting the 
result in (3.1), we may express the standard SMI weights W FT based on K 
fine-tuned snapshots as 


\V ft — 

= ii\C$caC h Y X s 

= /*(C")- 1 *5iC- 1 5 


(3.26) 


Let 4>ij l be the ij th element of Also since the desired signal is not 

present in this scenario, we are free to choose S = [0 1] T . When we make 
these substitutions in (3.26) and evaluate, we have 






c 0 

0 1 


4*12 


4* 21 


-1 

2 

1 


4*22 


C* 0 

0 1 


= M 
= ^ 


\c\ 2 4> n 

c 4* 12 

0 

c *<^21 

4*22 

i 


c 4*n 

4*22 


CWiCA 

W2CA 


(3.27) 


where w lCA , and w 2CA are the auxiliary and main element weights based 
on the corresponding course- adjusted (not fine-tuned) snapshots. In words, 
(3.23) and (3.27) state that the effect of fine-tuning the auxiliary element 
samples with the factor 1/c* is to change the resulting auxiliary element 


35 


weight by the factor c. Recall that this result assumes a particular scenario 
and S = [0 l] r . However, the sample scale factor 1/c* which depends on 
cable lengths and real component characteristics should be independent of 
signal scenario and choice of steering vector. 

Our goal is to determine the scale factor c. To do this, the single- 
interference scenario described above is set up and the uncalibrated (coarse- 
adjusted but not fine-tuned) system is used to estimate a coarse-adjusted 
auxiliary element weight Wica which does not yield a completely suppressed 
interference signal at the array output. By stepping the magnitude and 
phase of the auxiliary element weight through a range of values, we deter- 
mine the ideal fine-tuned weight Wipi which causes the output interference 
signal to be completely suppressed. We then calculate c = ft/wica (us- 
ing the result in (3.27)) and apply the fine-tuning factor 1/c* (see (3.23)) to 
the first auxiliary element samples in subsequent experiments. The analo- 
gous procedure is followed to fine tune the second auxiliary element samples. 

This calibration process may be fully implemented in software by us- 
ing the output interference power measurement capability described in the 
previous section as follows: 

Calibration Procedure 


I 


Set up the scenario 

A. 2-element array (main and first auxiliary). 

B. Only first interference signal, no desired signal, no added noise. 

C. Set attenuator a 4 of Figure 2.3 so that a significant 
amount of interference is injected into the main channel to be 
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cancelled with the interference in the first auxiliary. 


II Adapt the weight 

A. Take samples for snapshot formation (enough for steady state 
estimation). 

B. Estimate the first auxiliary element weight u;. Note: Since the 
weight is uncalibrated, it needs to be scaled in phase and magnitude. 

III Find phase and magnitude scale factors ( <f> 0 and M 0 ). 


A. 

B. 

C. 

D. 

E. 

F. 


Let M 0 — 1 and <p 0 — 0. 

Step the phase of the weight wM 0 e 3,>>0 through a range of values. 

At each step (varying </>): 

i. Apply the phase-adjusted weight e^wMoe 3 * 0 . 

ii. Calculate and save the interference power at the array output. 
Choose the phase (p that yielded the minimum output 
interference power, call it <p a . 

Set (j) 0 = <f>o + <Pa- 

Scale the magnitude of the weight t vM 0 e 3<t>0 through a 
range of values. At each step (varying M): 

i. Apply the magnitude-scaled weight MwM 0 e 3 ^° . 

ii. Calculate and save the interference power at the array output. 
Choose the magnitude scale factor M that yielded the minimum 
output interference power, call it M a . 


G. Set M 0 — M 0 M a . 

H. Repeat IIIB.-IIIG. until convergence of M 0 and <p a . 
IV In all subsequent experiments modify the first auxiliary 

element (scaled) samples by the factor e 3 *° /M 0 . 
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V Repeat I-IV using the second auxiliary element instead of the 

first in order to fine tune the second auxiliary element (scaled) samples. 


The above procedure has been made iterative because the VMOD’s do 
not give us completely independent control of the weight magnitude and 
phase. Let us now turn to the problem of generating a steering vector. 

3.7 Estimation of the Steering Vector 

The steering vector contains the information that allows the SMI array 
to distinguish the desired signal from the interfering signals. In order to 
estimate the steering vector for a given bench-generated signal scenario 
we begin by turning off the interference signals and noise generators so 
that only the desired signal and negligible system noise are present in the 
received signals. We then take snapshots and form the estimate 

1 K 

^ ~ x kl^k (3.28) 

n fc = i 

where X k is the k th snapshot and x fel is the first element of X k . We note 
in passing that (3.28) is simply the first column of the standard sample 
covariance matrix (3.4). Once the steering vector estimate is stored in 
memory, the interference signals and noise generators are restored to the 
desired levels. This method of steering vector generation will be used in 
later experiments. Throughout this report, however, the main antenna is of 
relatively large gain compared to the auxiliary elements. This assumption 
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allows us to approximate the steering vector as 


S = [0 0 1] T (3-29) 

where T denotes transpose and the third element of S corresponds to the 
main element. The subject of generating a steering vector when the re- 
ceived signals are from geosynchronous satellites is not addressed here. 
Now that the experimental modified SMI system has been described and 
all relevent procedures have been discussed, it is time to describe its per- 
formance through experimentation. 
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Chapter 4 

Experimental Results 


The experiments test the steady state performance of the modified SMI 
array system for various input signal scenarios. By steady state perfor- 
mance we mean that performance achieved by using many snapshots in the 
weight estimation. The performance is studied using a step-by-step ap- 
proach. First, the system is tested using standard SMI (i.e. with F — 0). 
These tests establish the invariance of the calibration factors to changes in 
signal scenario, compare the experimental results with theory, and point 
out a lack of suppression of weak interference. The modified SMI algo- 
rithm is then used to increase the suppression of weak interference while 
maintaining a strong desired signal. In all experiments, the noise power 
in each element of the array is (approximately) the same and the desired 
signal, when present, is incident from broadside and present only in the 
main channel. The procedure used for conducting an experiment is 

1. Set the attenuators of the array simulator so as to present certain 
desired and interference signal levels to the array processor. 
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2. Turn off the noise generators. 


3. Set the weights to zero. 

4. Estimate the desired and interference signal levels using (3.15)-(3.17) 
(for J = 1, 2, and 3) 

5. Set the noise generators to the desired levels. 

6. Estimate the noise power in each array element using (3.18) (for J = 
1, 2, and 3). 

7. Collect samples with which to form snapshots (making sure that 
enough samples are taken to yield near steady state weights). 

8. Estimate the modified SMI array weights. 

9. Apply the estimated weights. 

10. Estimate noise power in the array output using (3.18) with J=4. 

11. Turn off noise generators. 

12. Estimate the desired and interference signal powers in the array out- 
put (J=4) using (3.15)-(3.17). 

13. Calculate performance parameters such as SNR and IS1 of (3.19) and 

( 3 . 20 ). 
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4.1 Experiment Is Two elements, one inter- 
ference, standard SMI, vary interference 
power 

In the first experiment, each auxiliary channel is tested separately by using 
the two-element (main and a single auxiliary) capability of the experimental 
system and standard SMI ( F = 0). No desired signal is present. A single 
interference of varying power is incident on the array from a fixed direction. 
The direction of arrival (52° off broadside in this case) is unimportant since 
any element pattern considerations are absorbed into the received signal 
powers. INR(aux)-INR(main) is fixed at 5.5dB. Figure 4.1 shows the in- 
terference suppression versus INR(main) for experiment and theory. One 
of the experimental curves corresponds to the two-element array consisting 
of the main and first auxiliary elements whereas the other corresponds to 
the main and second auxiliary element combination. A good agreement 
between theory and experiment is observed. In the linear region of the 
plot, a 5dB increase in interference power leads to a lOdB increase in its 
suppression. In the region INR(main)=10~15dB the small output inter- 
ference power becomes difficult to measure even with the noise generators 
turned off, due to ever-present residual noise; the main sources of which 
are the detector amplifiers A5 and A6. Note that for weak interference 
(INR(main) <-5dB), standard SMI yields less than 7dB of interference sup- 
pression. 
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Interference 


auxl and maln,latt=5.5dB 
aux2 and main,latt=5.5dB 
theory for 2-element array, latt=5.5dB 


One interference, no desired, main and one aux. 



Input INR in main [dB],(Pn fixed) 


Figure 4.1: Experiment 1: Suppression of a single interference of varying 
power with a two-element array. 
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4.2 Experiment 2: Three elements, one in- 
terference, desired signal, standard SMI, 
vary interference arrival angle 

Rather than varying the power level of the single interference we vary the 
direction of arrival (DOA) in Experiment 2. This time, standard SMI and 
all three elements are used and SNR(main)~16dB. The desired signal level 
in the auxiliary elements is negligible. The purpose here is to test the phase 
shifters of the array simulator and the phase-shifting ability of the vector 
modulators (weights). In addition, the invariance of the calibration factors 
to different scenarios is tested. The fixed input interference signal levels 
appear in Table 4.2. Figure 4.2 shows the results. One set of curves (theory 
and experiment) corresponds to the first interference signal whereas the 
other pair corresponds to the second interference signal. Five data points 
were used to form each experimental curve. There is reasonable agreement 
between theory and experiment since only one data point is further than 
2dB from its theoretical value. As expected, interference suppression is 
theoretically and experimentally independent of the arrival angle (as long 
as desired and interference signals arrive from different directions). Also 
note that the level of interference suppression provided by the standard 
SMI algorithm for this particular scenario is about 7.5dB. This may not be 
enough in some applications. 
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Interference suppression [dB] 


int.#1 , theory 
int.#1 , exper. 
int.#2, theory 
int.#2, exper. 


Vary direction of arrival 



DOA (degrees from broadside) 


Figure 4.2: Experiment 2: Suppression of a single interference of varying 
arrival angle with a 3-element array and SNR(main)=16dB. 
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INR(main) 

INR(auxl) 


Int. 1 only 

-3.78 

1.21 

-16.84 

Int. 2 only 

-4.67 

-16.48 

1.47 


Table 4.1: Interference signal levels for Experiment 2. 

4.3 Experiment 3: Three elements, one in- 
terference, desired signal, standard SMI, 
vary interference power 

Experiment 3 is a more complicated version of Experiment 1. Standard 
SMI and all three elements are used to suppress the first interference 
signal whose power level is varied. The desired signal is present with 
SNR(main)=16.5dB. Again, the desired signal level in the auxiliary ele- 
ments is negligible. Figure 4.3 shows the results. For one pair of curves 
(theory and experiment), INR(auxl)-INR(main) is fixed at 6.5dB whereas 
for the other pair it is fixed at 9.7dB. The difference in interference sup- 
pression between the two curve pairs is the result of increasing the aux- 
iliary element gain. Thus, in the linear portion of the curves, a 3.2dB 
(=9.7-6. 5) increase in auxiliary element gain has led to approximately a 
5dB increase in interference suppression for a given main element sidelobe 
level. Note again the lack of interference suppression for weak interference 
levels and low auxiliary element gain. For example, from the lower set 
of curves in 4.3 (corresponding to the lower gain auxiliary elements) we 
see that the array provides less than 10 dB interference suppression when 
INR(main ) < —4 dB. Theory and experiment show good agreement in 
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Figure 4.3 especially in the middle range -8dB<INR(main)<8dB. Above 
this range it becomes difficult to measure the small output interference 
power. 

Figure 4.4 shows the output SNR curves for Experiment 3. Note that, 
as theory predicts, the SNR is maintained at the output. 

4.4 Experiment 4: Three elements, one in- 
terference, desired signal, standard SMI, 
vary main antenna sidelobe 

In the fourth experiment (Figure4.5) using standard SMI, the first inter- 
ference signal is present, the first auxiliary element gain (INR( auxl)) is 
fixed, INR(aux2)= INR(main)-12dB, and the main element sidelobe level 
is allowed to vary. The second interference is off and the desired signal 
is present with SNR(main)=16.06dB. The top pair of curves are theory 
and experiment for INR(auxl)=7.63dB, and the middle and lower pair are 
for INR(auxl)=4.20dB and 1.20dB, respectively. The discrepancy between 
theory and experiment in the top pair of curves is explained, again, by 
the difficulty in measuring such small output interference powers. In the 
measureable regions, however, results are quite good. It is interesting to 
note that when the interference level in the main element is small relative 
to that in the auxiliary element (i.e. the main element sidelobe level is 
small compared to the auxiliary element gain), the theoretical interference 
suppression is relatively independent of the sidelobe level. From Figure 

4.5 we see that for low-gain auxiliary antennas (INR(auxl)<1.20dB) and 
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Interference Suppression [dB] 


INRA1 - INRM 4- 6 5dB, exp. 

INRA1 - INRM + 9.7dB, exp. 

INRAI* INRM + 6.5dB, theory 

INRA1 = INRM + 9.7dB, theory 


3 ELEMENTS, II AND DESIRED, SNRM=16.50dB 



-15 -10 -5 0 5 10 15 

Input INR In main (dB), (Pn fixed) 


Figure 4.3: Experiment 3: Suppression of the first interference signal with 
a 3-element array, SNR(main)=16.5dB. The top pair of curves are theory 
and experiment for fixed INR(auxl)-INR(main)=9.7dB. The bottom pair 
are for INR(auxl)-INR(main)=6.5dB. 
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SNR [dB] 


SNRin 

SNRout.INRAI - INRM + 6.5dB,theory 
SNRin-4 

SNRout-4,INRA1 = INRM + 9.7dB, theory 

SNRout.INRAI = INRM + 6.5dB,exp 

SNRout-4,INRA1 = INRM + 9.7dB,exp 

3 ELEMENTS, II AND DESIRED, SNRM=16.50dB 



Input INR In main [dB], (Pn fixed) 


Figure 4.4: SNR plots for Experiment 3: The top three curves are in- 
put SNR(main), theoretical SNR(out), and experimental SNR(out) for 
INR(auxl)=INR(main) + 6.5dB. The bottom three curves are the same 
for INR(auxl)=INR(main) + 9.7dB and have been displaced -4dB for clar- 
ity. 
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weak interference in the main element (INR(main)<-5dB), there is less than 
7.5dB interference suppression. Figure 4.6 shows the input SNR(main) and 
the three SNR(out) curves for Experiment 4. There is little or no loss in 
SNR from input to output. This experiment and the last will help us to 
quantitatively assess the value of using modified SMI. 

4.5 Experiments 5 and 6: Three elements, 
one interference, desired signal, fixed in- 
put powers, modified SMI, vary fraction 
F 

In the above experiments, we observed a lack of weak interference suppres- 
sion in the case of low-gain auxiliary antennas. The purpose of Experiments 
5 and 6 is to observe interference suppression when modified SMI is applied 
to such cases. In both of these experiments all 3 elements are used and only 
the first interference and desired signal are present. The results are plots of 
interference suppression and SNR as the fraction F of (3.3) is varied. Table 
4.2 presents the input signal levels for the two experiments. INR(main) 
is fixed at -9.45dB in Experiment 5 and is -2.68dB in Experiment 6. Fig- 
ures 4.7 and 4.8 show the results of Experiment 5 and Figures 4.9 and 4.10 
are the results of Experiment 6. The SNR plots show input SNR(main) 
and theoretical and experimental SNR(out) and are presented “waterfall 
style” so that the curves do not lay atop one another (i.e. the SNR curves 
corresponding to the top pair of interference curves have been displaced 
-4dB). 
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Interference suppression [dB] 


i inral *7. 63dB, theory ilnra1-l.20dB, theory 

ilnra1-7.63dB, exper. i1nra1«l.20dB, exper. 

i1nra1=4.20dB, theory 
i1nra1=4.20dB, exper. 


One interference, desired SNRM=16.06dB, vary main sidelobe 

20 

18 
16 
14 
12 
10 
8 
6 
4 
2 
0 

-20 -15 -10 -5 



Input INR in main [dB].(Pn fixed) 


Figure 4.5: Experiment 4: Suppression of the first interference signal with 
a 3-element array, SNR(main)=16.06dB. The first auxiliary element gain is 
fixed and the main element sidelobe varies. The top pair of curves are theory 
and experiment for fixed INR(auxl)=7.63dB. The middle and bottom pairs 
are for INR(auxl)=4.20dB and 1.20dB, respectively. 
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SNR [dB] 


SNRin 

SNRout, i1nral-7.63dB, exper. 

SNRout, i1nra1=4.20dB, exper. 

SNRout, i1nra1=1.20dB, exper. 


One interference, desired SNRM= 16 . 06 dB, vary main sidelobe 
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Figure 4.6: SNR(main) and SNR(out) curves for Experiment 4. 
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INR(main) 


INR(auxl) 



Exp. 5 

-9.45 

16.50 

0.43 

-3.57 

Exp. 6 

-2.68 

16.50 

4.72 

0.53 


Table 4.2: Signal levels for Experiments 5 and 6. 

In all cases theory and experiment compare favorably. From Figures 
4.7 and 4.9 we see that interference suppression increases as the fraction 
F used in the modified SMI algorithm is increased from zero towards one. 
For example, in Figure 4.7 with INR(auxl)=0.43dB we see that interference 
suppression is increased 9.5dB by changing from F=0.0 (standard SMI) to 
F=0.8. One can also make the observation from Figure 4.7 that a modified 
SMI array with low-gain auxiliary antennas (INR(auxl)=-3.57) and F=0.6 
performs the same with respect to interference suppression as a standard 
SMI array with auxiliary antennas of 4dB (=0.43+3.57) higher gain. Sim- 
ilar observations can be made about Figure 4.9 of Experiment 6. From the 
SNR plots of Figures 4.8 and 4.10 we see that for F<0.9 there is less than a 
2.5dB loss of SNR from input to output. Experiments 5 and 6 have demon- 
strated that modified SMI yields increased suppression of weak interference 
with only a small loss in desired signal power. 
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Interference Suppression (dB) 
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Figure 4.7: Experiment 5: Suppression of the first interference signal versus 
fraction F using a 3-element modified SMI array with desired SNR(main) = 
16.50dB. INR(main)=-9.45dB. Bottom pair of curves has INR(auxl)= 
-3.57dB. Top pair has INR(auxl)=0.43dB. INR(aux2)=-20dB. 
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SNR (dB) 


SNRin 

SN Rout , IN R A1 =-3.57 .theory 

SNRout.INRAI =-3.57 

SNRin-4 

SNRout-4,INRA1 =0.43, theory 

SNRout-4,INRA1 =0.43 

Mod. SMI, INRM=-9.45dB, SNRM=16.50dB, f-gain tradeoff 

20 


15 


10 


5 



Fraction f 

Figure 4.8: SNR curves for Experiment 5. Top three curves are input 
SNR(main), theoretical SNR(out), experimental SNR(out) for INR(auxl) = 
-3.57dB. Similarly, the bottom three curves are for INR(auxl)= 0.43dB and 
have been displaced -4dB for clarity. 
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Interference Suppression (dB) 
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Figure 4.9: Experiment 6: Suppression of the first interference signal versus 
fraction F using a 3-element modified SMI array with desired SNR(main)= 
16.50dB. INR(main)=-2.68dB. Bottom pair of curves has INR(auxl) = 
0.53dB. Top pair has INR(auxl)=4.72dB. INR(aux2)=-17dB. 
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SNR (dB) 


SNRin 

SNRout.INRAI *0.53, theory 

SNRout.INRAI =0.53 

SNRin-4 

SNRout-4,INRA1 =4.72, theory 

SN Rout-4, INRA1 =4.72 

Mod. SMI, INRM=-2.68dB, SNRM=16.50dB, f-gain tradeoff 



Fraction f 

Figure 4.10: SNR curves for Experiment 6. Top three curves are input 
SNR(main), theoretical SNR(out), experimental SNR(out) to. = 'NR(«“1) 
0.53dB. Similarly, the bottom three curves are for INR(auxl) 
have been displaced -4dB for clarity. 
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Chapter 5 
Conclusions 


The problem of weak interference signals entering a satellite - ground sta- 
tion communication link through the sidelobes of the receive antenna has 
been studied experimentally. The receive antenna is made adaptive so 
that antenna pattern nulls can be formed in the directions of the interfer- 
ence signals. Standard adaptive antennas which are designed to maximize 
the output SINR fail to achieve adequate suppression of weak interference, 
especially when the auxiliary antenna elements have low gain. In a previ- 
ous solution to this problem, the feedback loops used to control the array 
weights were modified so that the two inputs to the feedback loop corre- 
lators had uncorrelated noise components. These modified feedback loops 
yielded weights which were less influenced by noise and as a result adapted 
to null the interference. This solution had the drawback of requiring either 
two antennas per auxiliary element, two amplifiers per auxiliary element, 
or both. 

In this report, an alternative solution using a modified SMI algorithm is 
studied experimentally. The modification involves subtracting a fraction of 
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the minimum eigenvalue of the estimated covariance matrix from its diag- 
onal elements before inversion. The modified covariance matrix resembles 
a standard covariance matrix formed in a less noisy environment and thus 
“frees” the weights to further suppress the interference. The SMI modifi- 
cation is made in software and thus has the advantage of not requiring any 
additional antennas or amplifiers. 

The modified SMI algorithm is implemented on a three-element experi- 
mental array system. The need for a weight (vector modulator) in the main 
channel is eliminated by normalizing the estimated weight vector so that 
the main channel weight is unity. Sampling of the main channel is accom- 
plished by sampling the array output while the auxiliary element weights 
are set to zero. The signal snapshots are the input to the SMI algorithm 
and are formed from scaled samples. Such scaling provides a coarse adjust- 
ment for the losses present in the experimental system. The modified SMI 
algorithm is presented in detail. In order to quantify the performance of 
the experimental system the signal powers in each channel are estimated 
from samples. This is made possible by a special pulse modulation scheme 
in which the desired and interference signals are separated in time. Finally, 
a fully-automated, iterative procedure for system calibration is presented. 
In this procedure, amplitude and phase scale factors for the auxiliary signal 
samples are determined by comparing the weights that null the interference 
in a simple signal scenario with those that are estimated by the algorithm. 

The experiments that are presented in this report focus on steady state 
performance (i.e. a large number of snapshots are used to form the weight 
estimate for any particular scenario). In all cases, the experimental results 
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compare favorably with theoretical predictions based on the ideal array 
model. As expected, the interference suppression increases as the power 
of the interference signal incident on the array increases. In particular, 
it is noted that interference suppression is relatively independent of the 
main antenna sidelobe level when that sidelobe level is small relative to 
the auxiliary element gain. In other words, the interference suppression is 
to a great extent determined by the gain of the auxiliary antennas in such 
situations. When low-gain auxiliary elements are used, standard SMI leads 
to only small amounts of interference suppression. When modified SMI is 
used the suppression is increased by as much as 13dB in the experiments. 
Thus, for large numbers of snapshots the experimental modified SMI array 
can provide the desired performance as predicted by theory. 
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Appendix A 

Experimental System Software 
(Fortran IV) 

C MASTER PROGRAM FOR RUHIRG ADAPTIVE ARRAY SYSTEM 

COHMOR/OUTDTA/X , AVE , AVCOR , OFF , 

I LHAG.HHAQ.LPH.HPH.ISC.CPHI 

COMMOI/CTST/AVCD , AVIF , AVISE , ST1 
C0MH0R/SCIR0/LATT,LPHI,P1 
COMMOI/WTBLK/LDAO , VI ,WQ , VMI , VMY , IW 

COHMOR/SPVRS/PVR , PISE , ITER , IFL ,OFLAG , I FLAG , AIS , AIST 
COHMOI/SMI/COV ,IC0V .ICOV1 .ICOVIO , ICT 

COMPLEX X(3 ,04) ,AVCOR(2) ,0FF(3) ,AVCD<2) ,AVIF<2) ,AVISE<2) 

COMPLEX COV ( 3 , 3) , CO VCOP (3 , 3 ) . COVIIV (3 , 3) , STEER (3) , ST1 
COMPLEX CW(3) ,EVEC(3,3) ,DET ,ID(3,3) ,STR,CVTBST,WEI0HT(2O,2) 

REAL VMX(40 ,4) , VMY (40 .4) .POWERS (20, 4) ,ICOV 
REAL EVAL(3) ,LVECT(3) ,MVECT(3) ,PIMII,EIGVAL(20,3) 

REAL PWR(3,3) ,P*SE(3) , PB(3) ,PA(3) ,AIS(3) , AIST,CPHI(2) 

REAL AVE(0,3) , LMAG ( 2 ) , HMAO ( 2 ) ,LPH(2) ,HPH(2) ,PI ,P(2) .SIZE 

REAL LPHI (2 , 2) ,PI(4) ,VI (2) ,WQ(2) .TEMPI (2) ,TEMPQ(2) ,WI0(2) ,WQ0(2) 

IVTEGER IDAC0,LATT(3) , LDAO , OFLAG , EV (4 ) ,ICT, ICF.IDIM 
IITEGER IC0V1, ICOVIO 

L0GICAL*1 FLIME(0) ,FTYPE(4) ,FLSPEC(1B) ,IHE(9,4) 

DATA FTYPE/ » . > , , »A» , *T 9 /CPHI/2*0 . 0/ 

DATA IDAC0/ M 170440/0FF/3* (0.0, 0.0)/ 

DATA VI/2*0.0/VQ/2*0.0/ 

DATA PI/3 . 141S927/IFLAG/0/ 

DATA ISC/l7/STEER/(0. 0,0.0) , (0.0, 0.0) , (1 .0,0.0)/ 

C Evading data for VHOD control 

CALL SCOPY( 1 VM1I . DAT 1 , HHE(1 , 1) ) 

CALL SCOPY( ’VM1Q.DAT * , IHE(1 ,2) ) 

CALL SC0PY(>VM2I.DAT\IME(1,3)) 

CALL SCOPY ( 9 VM2Q . DAT 1 , IME (1,4)) 

DO IB 1*1,4 

OPEI ( UIIT* 1 0 , IAHE=IME ( 1 , 1 ) , T YPE* ’ OLD \ FORM* > UIFORMATTED * ) 

DO 10 J*1 ,40 

10 READ (10 , EID=ll) VMX( J ,1) , VMY (J , I) 

11 IV(I)=J 
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CLOSE (UIIT*10) 

IB COITIIUE 

C Enable pulse generator and set THA*s to TRACK mode 
CALL IPOKE(IDACO, "40000) 

2B VX<l)-0.0 

VI(2)*0. 0 

wq(i)=o.o 

VQ(2)=0.0 
CALL JWT 

C Enter Array Simulator parameters 
CALL JWIHT 
TYPE * , * » 


TYPE ♦,* Main and AUX1 elements only Enter 1* 

TYPE * , * Main and AUX2 elements only Enter 2* 

TYPE *,' All three elements Enter 3* 


TYPE 97 B 
ACCEPT 921.ICF 
IDIH*2 

IF(ICF.BQ.3)IDIM=3 

C 

TYPE *,» > 

TYPE * , 'BITER IUHBER OF EIGEIVECTORS TO KEEP:* 

ACCEPT 921 , IEIG 
C 

TYPE * , * * 

TYPE 978 
ACCEPT 922.FRAC 
TYPE*,* * 

TYPE*,* TRAISIEHT RUI ? EiTER 7 IF YES, 0 IF 10:* 

ACCEPT 921, ICODE 
C 

TYPE * , * * 

TYPE *,* Enter 1 to output powers to PVRS.DAT else Enter 0 :* 
ACCEPT 921.IP0W 

IF< IPOW . EQ . 1 ) OPEV (UII T*10 , SAME* * PVRS . DAT » , TYPE* * IEV » , 

+ FORM* * FORMATTED } ) 

C 

40 TYPE 919 

PAUSE 

C Calculate beginning powers 
TYPE 879 
CALL JVISC 
DO 60 L»l,3 

60 PB(L)*PWR(L,3) 

PIIT-PWR ( 1 , 3)+PVR<2 ,3) 

TYPE 929 
PAUSE 
■FLAG*1 
TYPE 879 
CALL JVISC 
IFL-1 

CALL RDOTPT 

TYPE 928 

ACCEPT 921, IRC 

IF (IRC .EQ. 1) GOTO 40 
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TYPE 966 
ACCEPT 970.I3C 
TYPE 977 
ACCEPT 921, ICT 
TYPE 979 
ACCEPT 923, LIMIT 
DO 380 1*1,3 
DO 390 J*1 , 3 
C0V(I,J)=(0. 0,0.0) 

390 COITIIUE 

380 COITIIUE 

■COV1*0 

■covio*o 

ITER*1 

C Start adaptation loop 
TYPE 968 
TYPE 969 
TYPE 966 


376 IF (ICODE .IE. 1) GOTO 490 
TYPE 967 
IC0DE*O 

490 IF (ICODE. IE ,7)G0T0 491 
FRAC * 0.0 

C Store present weights and apply ** ro weights 

C so that the main channel may be sampled as the array output. 

491 DO 600 K*1 ,2 


TEMPI (K)*WI(K) 
TEMPQ(R)*9QU) 
WI(lO*0.0 
VQ (R)*0 .0 


600 COITIIUE 

CALL JWWT 

c Update the covariance matrix by 12*ISC snapshots 

C then reapply adapted weights. 

CALL COVAR(ISC ,CPHI , OFF ) 


type*, cov(l , l) , cot(1 , 2) , cov(l ,3) 
type*, cov(2 ,1 ) , cot ( 2 , 2) , cot ( 2 ,3) 
type*, cot(3,1) ,cov(3,2) , cov(3,3) 

DO 610 1*1 , IDIH 

IF (I .IE. 3)WI(I) “TEMPI (I) 

IF(I.IE.3)WQ(I)*TEMPQ(I) 

DO 820 J*1 , IDIM 

II*I 

JJ*J 

IF(ICF.Ep. 1 .AID. I, EQ. 2)11*3 
IFCICF.EQ.1 . AID . J . EQ . 2) JJ*3 
IF(ICF.EQ. 2)11*1+1 
IF(ICF.EQ.2)JJ*J+1 
COVCOP(I,J)*COV(II, JJ) 


620 COITIIUE 

610 COITIIUE 

CALL JWVT 

c Get eigvals and eigvecs of cov . 

CALL HEIGEI (COVCOP , EVEC , EVAL , IDIM) 
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typm *» i i*l ,idim) 

C IORMALIZE EIGEIVECT0R5 

DO 623 1=1, IDIH 
SIZE*0.0 
DO 624 J*1 , IDIN 
SIZE-SIZ£+CABS(EVEC(J,I))**2 
624 COITIIUE 

SIZE*SQRT (SIZE) 

DO 62G J=1,IDIH 
EVEC(J,I)*EVEC(J f I)/SIZE 
626 COITIIUE 

623 COITIIUE 

C IIVERT COVARIAICB MATRIX 

802 IF ( ICODE . EQ . 6) FRAC=FRAC+0 . 1 

IF ( ABS (1. O-FRAC) .LT. 0.05) GOTO 670 
C IF( ABS (1. O-FRAC) .LT. 0. 0B)FRAC=0 .0 

C IFCl.l-FRAC .LT. O.OB)GOTO 670 

804 DO 621 1=1 , IDIH 

DO 622 J=1 , IDIM 
C0VIIVd,J)*(0. 0,0.0) 

DO 626 K*1,IEIG 

COVIIVCI, J)«COYIIV(I, J)+EVEC(I,I) *COIJG(EVBC(J ,![))/ (EVAL(I)-FRAC 

+ *EVAL(IDIM)) 

626 COITIIUE 

622 COITIIUE 

621 COITIIUE 

C CALCULATE COMPLEI WEIGHTS CV 

DO 630 1*1, IDIM 

CV(I)*(0. 0,0.0) 

DO 631 J»1 , IDIH 

S TR“ STEER (J ) 

IFdCF.EQ.l .AID. J . EQ . 2) STR=STEER(3) 

IF( ICF . EQ . 2) STR*STEER( J+l ) 

CW(I)=CV(I)+COVIIY(I, J)*STR 

631 COITIIUE 

630 COITIIUE 

C IORMALIZE WEIGHTS TO CV(IDIM)=m«in ohuui«l vtight 

DO 632 1=1, IDIM 
CW ( I ) =CW ( I) /CW ( IDIM) 

632 COITIIUE 

C PICK OFF I AID Q COMPOIKITS OF WEIGHTS 

DO 601 K**l , 2 
TEMPI(K)=WI(K) 

TEMPQ(K)*WQ(K) 

601 COITIIUE 

602 DO 633 K=l,2 

WI(K)= REAL(CW(K) ) 

WQ(K)=-1 . 0*AIMAG(CW(K) ) 

633 COITIIUE 
IF(ICF,EQ.1)WI(2)*0.0 
IF(ICF .EQ . 1 )WQ(2)=0 .0 
IF(ICF .EQ . 2) WI (2)*WI(l) 

IF(ICF.EQ.2)WQ (2)=WQ(1) 

IF( ICF . EQ . 2 ) WI ( 1 ) =0 . 0 
IF ( ICF . EQ . 2 ) WQ ( 1 ) =0 . 0 


66 


DO 634 E=l,2 

IF (IBS (WI (R) ) . LE . 1 . 0 .AID. ABS (WQ(R) ) .LE . 1 . 0) GOTO 634 
TYPE 972, R,R 
VI(R)=TEMPI(R) 

WQ(R)=TEMPQ (R) 

634 CONTIIUE 

IF(ICODE. EQ.7)G0T0 637 
IF (MODCITER^) .IE. 0) GOTO 63S 

637 DO 636 R=l,2 

IF (K .EQ. 1) TYPE 884.ITER 
TYPE 885,1,1 

TYPE 887, TEMPI (R) ,TEMPQ(R) ,VI(R) ,WQ(R) 

TYPE * , • 1 
636 COITIIUE 

636 CALL JWVT 

IF(ICODE.Eq.6)GOTO 710 

IF(IPOW.EQ.O) GOTO 456 

TYPE 919 

PAUSE 

TYPE 879 

HFLAG=0 

CALL JWISC 

TYPE 929 

PAUSE 

IFLAG=1 

TYPE 879 

CALL JWISC 

IF(ICODE . EQ . 6) TYPE 891 , FRAC ,ICOV 

WRITE (10 , *)ICOV , FRAC 

WRITEClO,*) (PWR(I,3) ,1=1 ,3) ,PISE(3) 

WRITEUO,*) (CW(I) ,1-1 , IDIM) 

WRITE(10,*) (EVAL(I) ,I=1,IDIM) 

456 IF (IRSP .EQ. 0) GOTO 640 

IF (ITER .LE. 250) WRITE(IO) (WI(I) ,WQ(I) ,1=1 .2) , 

I (AVCOR(I) , 1=1 ,2) , (AVCD(I) ,1=1,2) , (AVIF(I) ,1=1 ,2) , (AVISE(I) ,1=1,2) 

I, ITER 

640 IMCH=ITTIIR() 

ICRLF=ITTIIR() 

ICRLF=ITTIIR( ) 

IF (IMCH .LT. 0) GOTO 800 
IF (IMCH .Eq. 83) GOTO 660 
TYPE 878 
GOTO 800 

650 IF(MOD( ITER ,5) . Eq . 0) GOTO 670 

TYPE 884, ITER 
DO 660 R=1 , 2 
TYPE 886, R,R 

TYPE 887 , TEMPI (R) ,TEMPQ (E) ,WI(R) ,WQ(R) 



TYPE 



660 

COITIIUE 



670 

TYPE 888 




TYPE 

Cent inue adaptation 

— Entar 1 * 


TYPE 

Exit adaptation loop- 

— Entar 2 9 


TYPE 

Calculate jaamar suppression--- 

— Enter 3* 


TYPE *, ’ 

Output intarmadiata raaulta 

Entar 4 5 
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TYPE Co input# calibration factors Enter 6 * 

TYPE Vary frac Enter 6* 

TYPE * , * lumbar of eigenvectors Entar 8 1 

TYPE *,* Assign staaring vector, use COV Enter 9* 

TYPE *,* Rasat staaring vector to 001 Entar 10' 


TYPE 880 

ACCEPT 881 , ICODE 
IF (ICODE .EQ. 1) GOTO 800 
IF (ICODE .EQ. 5) GOTO 801 
IF (ICODE .EQ. 6) FRAC=-0.1 
IF (ICODE .EQ. 8) GOTO 802 

IF (ICODE .EQ. 8) TYPE * , * EVTER IUMBER OF EIGEIVECTORS TO KEEP : ’ 
IF (ICODE .EQ. 8) ACCEPT 921, VEIG 
IF (ICODE .EQ. 8) GOTO 670 
IF (ICODE .EQ. 9)STBER(1 )=C0V(1 , 1) 

IF (ICODE .EQ. 9)STEER(2)=C0V(2 , 1) 

IF (ICODE .EQ. 9)STEER(3)=C0V (3,1) 

IF (ICODE .Eq. 9) GOTO 670 
IF (ICODE .EQ. 10)STEER(1)*(0. 0,0.0) 

IF (ICODE .EQ. 10)STEER(2)=(0. 0,0.0) 

IF (ICODE .EQ. 10) STEER(3)=(l .0,0.0) 

IF (ICODE .EQ. 10) GOTO 670 
IF (ICODE-3) 880,710,720 

710 TYPE 919 

PAUSE 
TYPE 879 
IFLAG=0 
CALL JWISC 
DO 71B 1=1,3 

715 AIS(I)=10.0*ALOG10(PB(I)/PWR(I ,3) ) 

AIST=10.0*AL0G10(PIIT/(PWR(1 , 3)+PWR(2 ,3) ) ) 

TYPE 929 

PAUSE 

IFLAG=1 

TYPE 879 

CALL JWISC 

TYPE 890, A I ST 

TYPE 889, (AIS(I),I=1,3) 

IF ( ICODE. EQ .6)TYPE 891 , FRAC .ICOV 

891 FORMAT ( 1 FRAC = \F8.S,' ICOV = J ,F7.0) 

IF ( ICODE. EQ.6)G0T0 802 
GOTO 670 

720 CALL RDOTPT 

GOTO 670 

800 IF (LIMIT ,LE .0)GOTO 803 
LIM=LIMIT 

IF(*C0V .GE. FLOAT (LIMIT) ) LIMIT=LIMIT+12*ISC 
IF(ICOV.GE. FLOAT(LIM) ) GOTO 670 

803 IF( ICODE . IE . 7) ITER=ITER+1 

IF (ICODE . IE . 7) GOTO 376 
IF ( FRAC . EQ . 0 . 9 ) ITER=ITER+1 
IF (FRAC. EQ. 0 . 9) GOTO 375 
IF(FRAC. EQ.0.7) FRAC=0.9 
IF(FRAC.EQ.O.O) FRAC=0.7 
GOTO 804 

801 5CBEST=1 . 0 
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DO 812 J*l,6 
PIMI 1*909999. 999 
DO 811 I»-10,10 
THET AD=FLOAT ( I ) 

THETA=PI/180 . *THETAD 
CI*COS (THETA) 

CQ*3II (THETA) 

CVTEST=CV(1)*CHPLI(CI,CQ)*SCBEST 
VI(ICF)* REAL(CVTEST) 

VQ(ICF)=-1 .0*AIMAG(CWTEST) 

CALL JVVT 
CALL JWISC 

IF ( PVR ( ICF , 3 ) . LT . PIMI V ) PHBEST*THETAD 
IF(PWR(ICF, 3) . LT.PIHII)PIMII*PWR(ICF ,3) 

TYPE* , THETAD , PHBEST 
811 COITIIUE 

TYPE 462 

462 FORMAT ( * EITER PHBEST BETVEEI -10. AID 10.0 

ACCEPT 924, PHBEST 
PIMI 1=9999999. 999 
5C=0 . 90 
DO 810 1=1,21 

IF( CABS( CW(1) *SC ) .GE. 1.0)G0T0 812 
CI=C0S(PHBEST*PI/180. ) 

CQ*SII(PHBEST*PI/180. ) 

CVTEST=CV ( 1 ) *SC*CMPLX ( Cl , CQ ) 


810 

463 

812 

860 


877 

878 

879 

880 
881 
684 


VI ( ICF ) = REAL (CWTEST ) 

VQ(ICF)=-1.0+AIHAG(CVTEST) 

CALL JVVT 
CALL JWISC 

IF(PWR(ICF, 3) . LT . PIMII ) SCBEST*SC 
IF (PVR (ICF, 3) . LT . PIMII )PIMII*PVR( ICF ,3) 

TYPE* , SC , 5CBEST 
SC=SC+ .01 
COITIIUE 
TYPE 463 

FORMAT ( 1 EITER SCBEST BETVEEI 0.9 AID 1.1 O 
ACCEPT 924, SCBEST 
COITIIUE 
GOTO 670 

IF (IRSP .EQ.l) CLOSE (UIIT=10) 

IF ( IPOW. EQ. 1) CLOSE (UIIT=10) 

TYPE 877 

TYPE *,’ End program execution Enta: 

TYPE 880 
ACCEPT 881 ,IECD 
IF (IECD .EQ. 2) GOTO 990 
IF (IECD .IE. 1) GOTO 860 
GOTO 26 

FORMAT (// 9 Conduct another test Ent 

FORMAT (/' IIADVERTAIT CHARACTER EITEREDV) 
FORMAT(/’ Poser levels being computed 
FORMAT ( J Enter code for desired option: > ,t) 
F0RMAT(I2) 

FORMAT (//10X , 1 Iteration number ’,13) 


69 


885 

887 

888 

889 

890 

918 

919 

920 

921 

922 

923 

924 

925 
928 

927 

928 

929 
933 
954 

936 

937 

964 

966 

966 

967 

968 

969 

970 

971 

972 

973 

974 
976 

976 

977 


978 

979 
990 
999 


FORMAT (/ 131, 'LAST tfC \I1 , » ) ’,17X, ’CURRENT V(\I1,’)>) 

FORMAT ( 61, 2 (El 0. 3,21) ,4X , 2CE10 . 3, 2X) ) 

FORMAT (/ / 1 0PTIOI5 : * ) 

FORMATC/’ Jamraer-i: »,F8.4,’ Jammer-2: ’,F8.4,’ Desired 
I Sig: \F8.4/) 

FORMAT (/ * IITERFEREICE SUPPRESSION *,F8.4/) 

FORMAT (/’ For phase calibration, set all leakage attenu 
lators to 70dB. ’) 

FORMATC/ ’ Remove noise and hit [RETURI] : ’ ) 

FORMATC/’ Enter code for storing weight values: ’,$) 

FORMAT (12) 

FORMAT (F3. 1) 

F0RHATCI6) 

FORMAT (F6. 3) 

FORMATC/’ Enter filename ( <*8 characters): ’,$) 

F0RHAT(Q,8A1) 

FORMATC//' To storeCon floppy) weight vector and its 
I progression’) 

FORMATC/’ Enter 1 to recalculate else enter 0: ',$) 

FORMATC/’ Restore noise and hit [RETURI] : 1 ) 

FORMATC/’ Remove interference and noise for steering vector’) 
FORMATC’ calculation and hit [RETURI]:’) 

FORMATC/’ Restore interference and noise and hit [RETURI]’) 
FORMATC/’ Enter 1 to compute and use steering vector 
X else enter 0: * ,|) 

FORMATC//’ Enter loop gain: ’,$) 

FORMATC//’ Enter number of scans to average over: ’,$) 
FORMATC//’ Beginning Adaptation...’/) 

FORMATC//’ Continuing Adaptation...’/) 

FORMATC//’ Type ”S" and hit [RETURI] to suspend iteration’) 
FORMATC’ at any point and go to OPTIOI menu’) 

F0RMATCI3) 

FORMATCFIO. 6) 

FORMATC/’ WEIGHT #’,I1,’ TOO LARGE FOR VM0D-’,I1,» TO HAIDLE’/} 
FORMATC’ AUX - 1 = ’,F7.2,3X,’ AUX-2* \F7.2/) 

FORMATC/’ Feedback loop differential phase shifts:’) 

FORMATC’ Enter array configuration code: ’,$) 

FORMATC/’ Steering vector being computed...’/) 

FORMATC/’ Do not include crossterms EITER 0’/ 

X ’ Include d-j and j-j crossterms EITER 1’/ 

X ’ Enter covariance calculation method: ’,$) 

FORMATC’ Enter value for FRAC : ’,$) 

FORMATC’ Enter # of snapshots to take C<=0 for no limit): ’,$) 
TYPE *, ’ THE BID’ 

STOP 

BID 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C Subprogram JWIIIT--set array simulator parameters 
C Program to test control of programmable attenuators of array 
C simulator using parallel output port(DRVll) 

SUBROUTINE JWIIIT 

COMMOI/SCIRO/LATT ,LPHI , PI 

IITEGER DRBUF , IVAL,LATT(3) , IM(3) , IP (3) 

REAL LPHI (2 , 2) ,PI(4) 
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DATA DRBUF/" 167772/ IBIO/ 15*0/ 


DATA IM,IP/3*0,3*0/X2/-l/ 

IVAL=0 

C Desired signal parameters: array index 1 

C Jammer 1 signal parameters: array index 2 
C Jammer 2 signal parameters: array index 3 

20 F0RHAT(//> Enter Main signal leakage attenuation(dB) : * p l) 

21 FORMAT ( * Enter Jainner- 1 , II , » leakage attenuation(dB) : » P I) 

25 FORMAT (13) 

30 FORMAT ( 1 Set Main signal leakage rotary attenuator to , f I2 p , dB») 

31 FORMAT ( * Set Jamner- 1 t II p 9 leakage rotary attenuator to 

32 FORMAT ( * and hit [RETURI]') 

40 FORMAT (' Enter phase shift from Jammer-' * p II , 

X * to MAII(degrees) : 1 ,$) 

41 FORMAT ( 1 Enter phase shift from Jammer-* ,11 . * to AUX- , I I1 P , : * P l) 

45 F0RMATCF8.2) 

50 FORMAT ( * Enter Main channel noise power(dBm): *,♦) 

51 FORMAT ( 1 Enter Correlator channel noise power(dBm): * ,1) 

52 FORMAT ( * Enter AUX-\I1,» signal channel noise pover(dBm) : * p |) 

55 FORMAT (F7. 2) 

DO 500 J*1 , 3 
K=J-1 

IF (J . EQ. 1) TYPE 20 
IF (J .IE. 1) TYPE 21 ,K 
ACCEPT 26,LATT(J) 

IH(J)*IIT(0. 1*LATT( J) -2 . 15)*10 
IF( IM( J) .LT. 0) IM ( J ) =0 
IP(J)*LATT(J)-IM(J) 

IF (J .EQ. 1) TYPE 30, IM(J ) 

IF (J .ME. 1) TYPE 31 ,K , IM( J) 

TYPE 32 
PAUSE 

IVAL=IVAL+IP( J)*2** (K*5) 

600 COITIIUE 

CALL IPORECDRBUF, IVAL) 

C Enter in leakage phase shifts, noise povers that are manually set: 

DO 600 J-l , 2 
R2=-l .0*K2 
TYPE 40 P J 

ACCEPT 46, LPHI(J ,1) 

TYPE 41 , J, J+K2 
ACCEPT 45,LPHI(J ,2) 

600 COITIIUE 

TYPE *,* 9 

DO 700 J-l ,4 
IF (J .EQ. 1) TYPE 50 
IF (J .EQ. 4) TYPE 61 

IF (J .IE. 1 .AID. J .IE. 4) TYPE 52, J-l 
ACCEPT 55 , PI( J) 

700 COITIIUE 

RETURI 
EID 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
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C SUBPROGRAM TO SCAI AID A/D COIVERT SAMPLES FROM A WAVEFORM 

SUBROUTIIE JWSCAI 
COMMOI/OUTDTA/X, AVE 
COKMOI/WTBLR/LDAO 
COMPLEX 1(3,64) 

REAL D(6,64) ,AVE(6,3) 

IITEGER ADCSR, CLCSR , ADDBR , CLBPR , KCSR , I AD ( 6 , 64 ) , IDACO , LDAO 
DATA ADCSR f ADDBR , CLCSR , CLBPR/" 170400 , M 170402 , "170420, "170422/ 

DATA IDACO/ "170440/ 

DO 10 IFX=1 ,6 
DO 9 JFI=1 ,3 

9 AVEClFI , JFX) =0 . 0 

10 COITIIUE 

11 FORMAT (IB) 

C Configure clock for external start (5T2) and single interval mods 
KCSR=8200 
ISTO»"40000+LDAO 
ICIT2=~63 

CALL IPOKE(CLBPR, ICVT2) 

C Using TTL bit from DACO holding ragistar to start pulsa ganarator (ACT. LOW) 

C and put THA's to track mods 
CALL IPOKE( IDACO, ISTO) 

CALL IPOKE(CLCSR.KCSR) 

C Loop to chock A/D dona bit and road converted valuta 
DO 100 1*1,64 
50 IC* IPEEKB (CLCSR) 

C Chock clock overflow bit, than flag ovarrun bit 
IF (IC .LT. 128) GOTO 50 
IF (IPEEKB ("170421) .AID. "20) GOTO 900 
CALL IPOKE( CLCSR, KCSR) 

CALL I POKE (CLBPR , ICIT2-J) 

DO 95 K*1 , 6 

IASR*R*256+1 

CALL IPOKE(ADCSR, IASR) 

C Chock for A/D dona bit bafora reading converted data 
60 IF (IPEEKB (ADCSR) .LT. 128) GOTO 60 

C IF (IPEEKB ("170401) .GE. 128) GOTO 910 

IAD(R,J)=IPEEK (ADDBR) 

96 COITIIUE 

C Reset THA's back to TRACK 

CALL IPOKE( IDACO, LDAO) 

CALL IPOKE(IDACO.ISTO) 

100 COITIIUE 

C Disable clock starts, leaving pulsa gan running 
CALL IP0KE(CLCSR,0) 

CALL IPOKE( IDACO, LDAO) 

C Convert data from offset binary 
DO 200 12*1 , 64 
DO 195 R2*l ,6 

C SCALE ACCOUITS FOR LOSSES THROUGH VMDDS .POWER DIVIDERS 
C AID DIFFEREIT AMPLIFIER SETTIIGS AT THE A/D COIVERTERS. 

IF(K2.LE. 4)SCALE=0 . 36481 
IF(K2 . GT. 4)SCALE=1 .4142 
D(K2, J2)=(IAD(K2, J2)-"4000)*0.0025*SCALE 
C TYPE ♦,K2,J2,D(K2,J2) 
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10B COmiUE 

C— CALIBRATE (FIIE-TUIE) THE SAMPLES 
ALPHA*D(1,J2) 

BETA =D(2.J2) 

D(1,J2)=(- . 8761 8)* ALPHA- ( . 33075) *BET A 
D(2 , J2)=( .33075) * ALPHA +(-. 87618) +BETA 
ALPHA*D(3, J2) 

BETA -D(4 , J2) 

D(3 , J2)- (" . 43634) * ALPHA- ( . 87702) *BETA 
D(4»J2)=(. 87702) * ALPHA +(-. 43634) *BETA 

cnv=-i.o 

DO 199 J*l,3 
IF (J .EQ. 3) CIIV*1.0 

X ( J , J2)=CIIV*CMPLX(D(2*J-1 , J2) ,D(2*J,J2)) 

199 COITIIUE 

200 COITIHUE 

C Averaging pulse levels for comparison 
C MAII--from samples 2-13 
DO 4 BO 1*1,6 
DO 300 J2=2 ,13 

300 AVE(I ,3)=AVE(I ,3)+D(I,J2) 

C Jl--from samples 17,28 
DO 360 J2*17 ,28 

3B0 AVE(I , 1 )=AVB( I ,l)+D(I,J2) 

C J2 — from samples 33,44 

DO 400 J2=33,44 

400 AVE(I , 2)=AVECl ,2)+D(I,J2) 

AVE(I,l)=AVE(I,l)/12. 

AVE(I,2)=AVE(I,2)/12. 

AVE(I,3)=AVE(I ,3)/12 . 

450 COITIIUE 

GOTO 999 

900 TYPE FLAG OVERRUI BIT SET II KWCSR * , * J ® * » J 
999 RETURI 

EID 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c Program to calculate signal and noise power at array output, 
C from samples taken with JVSCAI, then calculate J/I ratios 
C and the interference suppresion 
SUBROUTIIE JV1SC 
COHHOI/OUTDTA/X , AVE , AVCOR , OFF 
COMHOI /SPURS /PVR , PISE , ITER, IFL , OFLAG , I FLAG 
COMPLEX 1(3,64) ,0FF(3) ,Y ,AVC0R(2) 

REAL PWR(3,3) ,PISE(3) , AVE (6, 3) ,S (3,3 ,2) ,0FFR(3) ,0FFI(3) 
IITEGER IFL, OFLAG, IFLAG 
C Enter number of samples 

20 FORMAT(/ * Enter number of scans to average over 1 ) 

21 FORMAT (’ for power calculations: *,$) 

25 F0RMAT(I3) 

TYPE 20 
TYPE 21 
ACCEPT 25 ,IUM 
C MUM=100 

ISMP=5UM+12 . 0 
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I5HP2=IUM*18 .0 
IF (IFLAG .EQ. 1) GOTO 863 
DO 16 J=1 ,3 
DO 10 K=l,3 
PWR( J ,K)=0 .0 
S(J,K,1)=0.0 
S(J,K,2)=0.0 


10 COITIIUE 

0FF(J)=(O. 0,0.0) 

16 COITIIUE 

DO 810 J=1,IUH 
CALL JVSCAI 
DO 800 KCH=1,3 
DO 90 K=47,82 

90 OFF (KCH) =OFF(KCH)+X (KCH , K) 

60S DO 850 111=1 , 3 

DO 800 11=1,2 
L=2*RCH-2+H 

5(I!1,KCH,I1)=S(IM,KCH,M)+AVE(L, IM) 

800 COITIIUE 

860 COITIIUE 

800 COITIIUE 

810 COITIIUE 

DO 860 H=1 , 3 
0FF(M)=0FF(M)/ISHP2 
OFFR(M)=REAL(OFF(M) ) 

OFFI (M) =AIHAG (OFF 00) 

DO 840 1=1,3 
S(I,H,1)*S(I,H,1)/IUM 
S(I,H,2)=S(I t N,2) /IUM 

PWR(I,H)= (S(I,M, 1 ) -OFFR(M) )**2+(S(I,H, 2) -OFFI (H) ) **2 
840 COITIIUE 

850 COITIIUE 

GOTO 900 

C Boil* power calculations over portion of PRP whan no signal present. 

863 DO 864 M=l,3 

864 PISE(M)=0. 0 

666 DO 890 1=1, IUM 

CALL JVSCAI 
DO 888 H«i,3 
DO 870 K*47 ,82 
Y=I(H, K) -OFF (M) 

PISE(M)«PISE00+CABS(Y)**2 


870 COITIIUE 

888 COITIIUE 

890 COITIIUE 

891 DO 896 H=l,3 

896 PISE(M)*PISE(M) /ISHP2 *3.0 

C THE FACTOR OF 3.0 ABOVE TAKES DUTY CYCLE IITO ACCOUIT 

900 TYPE 910, (PVR (1,1) ,1*1,3) 

TYPE 911, (PWR(I ,2) , 1=1 ,3) 

TYPE 912, (PVR (I ,3) ,1=1 ,3) 

TYPE 914, (OFF(M) ,M=1 ,3) 

IF (HFLAG .EQ. 1) TYPE 913, (PSSE(I) ,N=1 ,3) 

910 F0RHATC/6I, >AUX-1 POWERS; * ,3(2X, E12 . 6) ) 
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911 FORMAT(/BX, f AUX-2 POWERS: » ,3(2X,E12 . 5) ) 

912 FORMAT(/BX, 'ARRAY OUT POWERS ,3(2X , E12 . 6 ) ) 

913 FORMAT(/BX, 'BOISE POWERS: \ 3(21, E12 . B) ) 

914 FORMAT (/BX, ’OFFSETS : * ,6E10.3) 

HFLAG“0 

9999 RETURB 

BSD 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

SUBROUTISE COVAR(BSC,CPHI f OFF) 

COHMOH/SMI/COV ,BCOV ,HCOVl ,HC0V10 ,ICT 
COHMOH/OUTDTA/X , AVE 
COMMOS /WTBLK/LDAO 

IBTEGER SCO VI , SCO VI 0 ,LDA0 , BSC , I , J , K , ICT 
REAL Cl , CQ , CPHI ( 2 ) , AVE ( 6 , 3) , RV , PI , ICOV 

COMPLEX COV (3, 3) , SS (3) ,1(3,64) ,OFF(3) 

COMPLEX PHRAH1 ,PHRAS2 ,PHRAS3 
DATA IRAB1/0/IRAS2/0/ 

PI=4*ATAS(1 . 0) 

c CHARGE COV FROM A COVARIABCE MATRIX TO A SUM OF 

C SBAPSHOT OUTER PRODUCTS 

DO 8 J=1 , 3 
DO 9 K=1 ,3 

COV( J , K)*COV ( J ,K) ♦( FLOAT (SC0V1 )+FLO AT (ICOVIO) *10000. ) 

9 COBTISUE 

8 COBTISUE 

IF( ICT . EQ . 0) JMAX=44 
IF(ICT.EQ.1)JMAX=12 
CI=1.0 
CQ-0 , 0 

C THE 100 LOOP ISCREASES THE IUMBER OF SIAPSHOTS OVER WHICH THE 

C COVARIABCE MATRIX COV IS AVERAGED BY ISC*12. 

DO 100 1=1, BSC ! SCAB BUMBER 

CALL JWSCAB !GET SAMPLES ACROSS ORE PRP, STORE IB X ARRAY 

C the 10 LOOP ADDS 12 SBAPSHOT OUTER PRODUCTS TO COV. 

DO 10 J=»l , JMAX ! * SBAPSHOT ' BUMBER FOR SCAB I 
C IEED 3 PHASE-AT-SAMPLE-TIME RV'S FOR EACH SBAPSHOT 

C 


C 

C 

C 

C 

C 

c 


(OBE FOR EACH SIOBAL COMPOHEIT) 

RV=2*PI*RAB(IRAH1 ,IRAB2) 

PHRAB1 “CMPLX ( COS ( RV ) , S IB ( RV ) ) 

IF( ICT .BQ .0) GOTO 11 
RV= 2 *PI*RAB (IRAB1 ,IRAB2) 

PHRAH2=CMPLX(C0S(RV ) , SIB (RV) ) 

RV=2*PI*RAB(IRAB1 , IRAB2) 

PHRAB3*CMPLX(C0S (RV) , SIB (RV) ) 

•THE 20 LOOP CREATES A SBAPSHOT VECTOR 
WHICH LEADS TO CROSSTERMS 
DO 20 Ml .3 'SBAPSHOT ELEMERT BUMBER 
-ADD IB THE DESIRED, Jl, ABD J2 COMPOBEBTS OF THE 
SIGHAL SBAPSHOT VECTOR SS . 

-GIVE EACH SIGHAL COMPOHEST A RAHDOM PHASE AT SAMPLE TIME 
WHILE MAIBTAIBIHG IBTERELEMEBT PHASE RELATIOBSHIPS . 

SS(JO = (0. 0,0.0) 

55 (H) =SS (K)+(X(R , 1+J) -OFF(K) ) *CMPLX (Cl ,CQ) *PHRAB1 
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SS(K)*SS (K)+(X(K f 10+J)-OFF(K) ) *CHPLX(CI ,CQ) *PHRAI2 
SS(K)=5S(K)+(I(K , 32+J) -OFF(K) ) *CHPLX(CI ,CQ)*PHRAI3 
20 COITIIUE 

IF(ICT.EQ.1)G0T0 12 

11 IF(J.EQ.1)G0T0 10 

IF(13.LT. J .AID. J.LT. 17) GOTO 10 
IF(28.LT.J .AID. J.LT.33)G0T0 10 

C THE 26 LOOP CREATES A SIAPSHOT VECTOR 

C WHICH LEADS TO 10 S|Ji OR J1|J2 CROSSTERMS 

DO 26 K=l,3 ! SIAPSHOT ELEKEIT IUMBER 
SS (K) = ( X(K f J)-OFF(R) ) +CMPLX (Cl ,CQ) +PBRAI1 
26 COITIIUE 

c LOOPS 30 AID 40 ADD A SIIGLE SIAPSHOT OUTER PRODUCT 

C TO COV. 

12 DO 30 L-l ,3 

DO 40 M=l,3 

COV (L ,M)-COV(L,M) + SS (L)*COIJG(SS (M) ) 

40 COITIIUE 

30 COITIIUE 

10 COITIIUE 

100 COITIIUE 

c UPDATE THE IUMBER OF SIAPSHOTS USED II THE COV MATRIX AVERAGE 

■C0V1 * IC0V1 + 12+ISC 

IF(IC0V1 .GE . ioooo)icovio=icovio+i 

IFdCOVl .GE . 10000) IC0V1=IC0 VI -10000 


IC0V=FL0AT(IC0V1 ) +FLOAT(ICOV10) * 10000 . 
TYPE *,ICOV,ISC 


c DIVIDE THROUGH BY THE IUMBER OF SIAPSHOTS II THE AVERAGE. 

DO 60 J=1 ,3 
DO 00 R=l,3 

COV( J , K) =COV ( J ,R) / (FLOAT (ICO VI ) + FLOAT(ICOVIO) *10000. ) 

00 COITIIUE 

60 COITIIUE 

RETURI 
EID 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTIIE HEIGEI (H,EVEC, EVAL.I) 

COMPLEX H(3 f 3) f EVEC(3,3) ,B(6,6) ,V(6) 

COMPLEX PUII (6,0) ,ALAM,P1 (0,0) , A(6,6) 

REAL ID(0,6) ,1,0(36) ,R(36) ,EVAL(3) , ALPHA 
DATA B/36* (0 . ,0 . ) / , W/0+ (0. ,0 . )/ 

DATA PUII/30* (0 . , 0 . )/, ID/30 *0 . / 

DATA Pl/30* (0. ,0. )/ ,C/30*0 . / 

DO 10 1*1,1 
ID(I ,I)=1 . 

DO 10 J*1,I 
PI ( I , J )=ID(I , J ) 

10 B(I,J)=H(I,J) 

DO 20 J=1,I-1 
J5=0 

DO 24 I=J+1,I 

24 IF(B (I , J) . IE . (0. ,0.))J5=1 

IF(J6 . EQ . 0)GOTO 20 
DO 16 1=1,1 
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IB 


30 
1 BO 

BO 

60 

40 

100 

80 

110 

90 

20 

70 

130 

120 


DO IB Kl-1,1 
PUII(I,K1)*ID(I,K1) 

K-0 . 

DO 30 I*J+1,I 
K=CABS (B(I , J) ) **2+K 
K=SQRT(K) 

JB=0 

DO 1B0 I=J+2,» 

IF(B(I,J) .HE. (0. ,0 . ) ) J6=l 

IF((JB.EQ.O) . AMD . ABS(K-REAL(B(J+1 , J) ) ) «LT . 0.1E-10) GOTO 20 
ALPHA=2 . *R**2-2 . *K*REAL(B( J+l , J) ) 

ALAH=ALPH A/ ( K* * 2 -K*B ( J+l , J ) ) 

V(J+1)=B(J+1,J)-K 

DO BO I=J+2 ,H 
V(I)=B(I,J) 

DO 60 I=J+1,I 
DO 60 K1=J+1 ,1 

PUII(I,Rl)=ID(I,Kl)-(ALAM/ALPHA)*tf(I)*COIJG<W(Ri)) 

DO 40 1=1 ,H 
DO 40 Kl=l , B 
A(I,K1)=(0. ,0.) 

DO 100 1=1,1 
DO 100 L=1,R 
DO 100 M=i ,1 
DO 100 111=1,1 

A(I,H1)=PUII(I,L)*B(L,H)*C0IJG(PU»I(M1,H))+A(I,H1) 

DO 80 1=1,1 
DO 80 Xl-1,1 
B(I ,K1 )=A(I , Rl) 

A(I,K1)*(0. ,0.) 

DO 110 1=1 , H 
DO 110 M1=1,I 
DO 110 M2=l ,1 

A (I ,M2)=P1(I ,M1) *COBJG(PUII (M2 ,H1 ) )+A ( I,H2) 

DO 90 1=1, S 
DO 90 111=1,1 
P1(I,M1)=A(I,H1) 

COITISUE 

Kl=l 

DO 70 1=1,1 
DO 70 J=1,I 
C(K1 )=REAL(B( J , I) ) 

K1=K1+1 

CALL EIGEI(C,R,»,0) 

Kl=l 

DO 130 1=1,1 
DO 130 J=1,I 

EVEC (I , J) =CHPLX (0.00,0.00) 

PURI (J ,I)=R(K1) 

Kl=Rl+i 
DO 120 1=1,1 
DO 120 J*1,B 
DO 120 Ml =1,11 

evec(i,mi)=pi(i,J)^pusi(J,mi)+evec(i,mi) 

DO 140 1=1, ■ 
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Rl=I+(I*I-I)/2 
140 EVAL (I)=C(K1 ) 

RETUR1 

EID 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

SUBROUTIIE EIGEV( A,R,V ,MV) 

REAL A (30) , R(30) , RABGE , FH , ABORM, AIORMX ,THR, I, Y 
REAL SIBX ,SIBX2 ,0051,005X2, SIICS 
5 RABGE»l.OE-0 

IF(HV-1 ) 10,26,10 
10 IQ=~B 

DO 20 J=1 ,1 
IQ-IQ+B 
DO 20 1=1,1 
IJ=IQ+I 
R(I J )=0 . 0 
IF(I-J)20,16 ,20 
16 R(IJ)=1 . 0 

20 COBTIBUE 

26 ABORM=0.0 

DO 36 1=1,1 
DO 36 J=I,B 
IF(I-J)30 ,35 ,30 
30 IA=I+(J*J-J)/2 

AB0RM=AB0RM+A(IA) *A (IA) 

36 COST HUE 

IF( AflORM ) 106,106,40 
40 AB0RM*1 . 414*SQRT ( ABORM) 

FB=1 .0*1 

ABRMX=ABORM*RABGE/FB 

I*D=0 

THR= ABORH 
46 THR=THR/FI 

60 L=1 

66 M=L+1 

00 MQ=(M*H-H)/2 

LQ=(L*L-L)/2 
LH=L+HQ 

02 IF(ABS (A(LH) )-THR >130,06,06 

06 IBD=1 

LL=L+LQ 

MM=M+MQ 

X=.6*(A(LL)-A(MM)) 

08 Y=-1*A(LM)/SQRT(A(L«)*A(LH)+X*X) 

IF(X)70,76,76 
70 Y=-Y 

76 SIIX=Y/SQRT (2 . 0* (1 .0+(SQRT(l . 0-Y*Y) ) ) ) 

SIBX2=SIBX*SIBX 
78 COSX=SQRT (1 .0-SIBI2) 

C0SI2=C0SX*C0SX 

SISCS=5IBX*C0SX 

ILQ=B*(L-1) 

IMQ=H* (H-l ) 
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DO 126 1=1,1 
IQ=(I+I-I)/2 
IF(I-L)80,116,80 
80 IF(I-H) 86,115,90 

85 IM=I+MQ 

GO TO 95 
90 IH=M+IQ 

95 IF ( I -L) 100, 106 , 105 

100 IL=I+LQ 

GO TO 110 
106 IL=L+IQ 

110 X=A ( IL)*COSI - A (IH) *SIII 

A (IM)=A(IL) *SIMX+A(IM) *COSX 
A(IL)=X 

116 IF(MV-1 ) 120, 125,120 

120 ILR=ILq+I 

IMR=IMQ+I 

X=R(ILR)*COSX-R(IMR)*SIIX 
R(IMR) =R(ILR)*SIIX+R(IMR) *COSX 
R(ILR) =X 
125 COSTIRUE 

X=2.0*A(LH)*SIICS 
Y=A(LL)*C0SX2+A(MM)*SIHX2-X 
X*A(LL)*SIIX2+A(MM) *C0SX2+X 

A(LH)=(A(LL) -A(MM)> *SIICS+A(LK>* (C0SX2-SIIX2) 

A(LL)=Y 

A(MM)=X 

130 IF(H-I) 136,140,135 

136 H=H+1 

GO TO 90 

140 IF(L-CI-I) >146,160, 146 

145 L=L+1 

GO TO 66 

1 60 IF( IHD-1 >160,165, 160 

166 IID=0 

GO TO 50 

160 IF(THR-AIRMI) 166 ,165,46 

165 IQ=-I 

DO 185 1=1, R 

iq=iq+R 

LL=I+(I*I*-I)/2 

jq=R*(I-2) 

DO 186 J=I,R 

jq=jq+i 

MM=J+(J*J-J)/2 
IF(ACLL)-A(KM>) 170,185,186 
170 X=A (LL) 

A(LL)=A(MM) 

A(MH)=X 

IF(MV-1)176, 186 ,176 
176 DO 180 X=1,H 

ILR=IQ+K 
IMR=jq+R 
X=R(ILR> 

RCILR)=R(IMR) 

180 R(IMR)=X 
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185 COITIIUE 
RETURI 
EID 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C Subroutine to print /type data from a particular experiment (JVMSV2) 
SUBROUTIIE RDOTPT 
COMMOI/OUTDTA/X , A VE , AVCOR , OFF 
X , LM AG , HMAG , LPH , HPH , ISC 
C0MH0I/5CIR0/LATT ,LPHI , Pi 

COMHOI/SPWRS/PVR , PISE p ITER , IFL , OFLAG , I FLAG , AIS , AIST 
COMPLEX 1(3,64) f AVC0R(2) # 0FF(2) 

REAL AVEC6.3) ,LMAG(2) ,HMA0(2) ,LPH(2) ,HPH(2) ,PI , AIS(3) , AIST 

IITEGER LATT (3) , OFLAG , IFLAG 

REAL LPHI (2 , 2) , PI (4) , PWR(3 f 3) , PISE(3) 

KD=-1 

ICR=13 

ILF=10 


15 


20 

25 


34 

35 

40 

41 


100 

101 

105 

106 

108 

109 


TYPE 111 

TYPE * i ’ Select Output Media: 1 

TYPE *,’ Monitor display 

TYPE 20 

ACCEPT 41, OFLAG 


FORMAT ( 

TYPE *, 

TYPE * , ’ Select Output : * 

TYPE * , 1 Signal Scenario-Simulator parameters — Enter 1* 
TYPE Signal /Jaimner and Signal/Ioise ratios 1 

TYPE * , ’ at each detector Enter 2* 

IF (IFL .EQ. 1) GOTO 34 

TYPE Sampled data from last scan Enter 3* 

TYPE DC Offsets, Correlation range for* 

TYPE last iterat ion(ave over ISC scans) Enter 4* 

TYPE * , ’ Interference suppression Enter 5 * 

TYPE * , ’ Change output media Enter 6* 

TYPE Done-Return to Main Program Enter 7’ 

FORMAT ( ’ Enter code for desired output: ’,$) 

FORMAT (II) 

TYPE 40 


Enter output media code: 


ACCEPT 41 , IPCODE 
IF (IPCODE .BQ. T) GOTO 9999 
IF (IPCODE .EQ. 6) GOTO 15 
IF (IPCODE .EQ. 5) GOTO 700 
IF (IPCODE .EQ. 4) GOTO BOO 
IF (IPCODE- 2) 150,970,950 
TYPE 111 

FORMAT (’ Main leakage attenuation= ’ ,12,’ dB’,2Al) 

FORMAT ( ’ Jammer- ' ,11 , 9 leakage atten.= ’,I2,» dB ’,2A1) 

FORMAT ( ’ J\I1,’/MAII leakage phase shif t- * ,F7 . 2 , ’ degrees’ 

X , 2X , 2A1 ) 

FORMAT (’ J ’ ,11 , ’/Aux-’ ,11 , ’ leakage phase shift* ’ ,F7. 2 , ’ degrees’ 
X.1X.2A1) 

FORMAT ( ’ Main channel noise power* 1 ,F6 . 2 , ’ dBm’,2X,2Al) 

FORMAT(’ Corr. channel noise power* 1 ,F6 . 2 , ’ dBm’,lX,2Al) 
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F6.2, 3 dBm 3 ,2Al) 


110 FORMAT ( 3 lux- 3 ,11 , 3 channel noise power* 3 , 

111 FORMAT*/) 

1BO IF (OFLAG .IE. 1) GOTO IB 

TYPE 100,LATT(1) 

DO 190 J*2 , 3 

190 TYPE 101 , J-l ,LATT(J) 

TYPE 111 
DO 220 J=l,2 
KD= (-1 ) ♦*( J~1 ) 

TYPE 106,J,LPHI(J,1) 

TYPE 106 ,J , J+RD,LPHI(J ,2) 

220 COITIIUE 

TYPE 111 
TYPE 108,PI(l) 

TYPE 110,1 ,PI(2) 

TYPE 110,2, PI(3) 

TYPE 109,PI(4) 

GOTO 26 

C Program section to output dc off Bat, correlation range, etc. 

600 IF(OFLAG .IE. 1) GOTO IB 

TYPE 111 
TYPE 910, ISC 
TYPE 111 
DO BSO J-l , 2 

BBO TYPE 906, J ,OFF(J) 

TYPE 911 , 0FF(3) 

DO 600 J=1 , 2 
TYPE BOB , J , AVCOR( J ) 

TYPE 901 , J ,LPH CJ ) 

TYPE 902 , J ,HPH(J ) 

TYPE 903 ,J ,LMAG( J) 

TYPE 904 ,J ,HMAG( J) 

TYPE 111 

600 COITIIUE 
GOTO 26 


C Output intarfaranca calculations 


700 


IFCOFLAG .IE. 1) GOTO 16 

TYPE 111 

TYPE 878 

TYPE 111 

TYPE 880, AIST 

TYPE 882,(AIS(J),J*1,3) 

TYPE 111 


GOTO 26 


878 

FORMAT 

( 3 

880 

FORMAT 

( 3 

882 

FORMAT 

( 3 

X 

3 Desired 

901 

FORMAT 

< 3 

902 

FORMAT 

( 3 

903 

FORMAT 

( 3 

904 

FORMAT 

( 3 

906 

FORMAT 

( 3 

906 

FORMAT 

(' 

910 

FORMAT 

( 3 


After 3 ,13,* iterations: 
IITERFEREBCE SUPPRESSIOI: 
Jammer-1: dB 3 , 3 


\2A1) 

3 ,F8.4,2A1) 

Jammer -2 : \F8.4, 3 


dB 3 


SMALLEST AUX- 3 , II, 3 CORR. PHASE* 3 ,F8.3, 2A1) 
LARGEST AUX- 3 , II , 3 CORR. PHASE* 3 ,F8 . 3, 2A1) 

SMALLEST AUX“ 3 , II , 3 CORR. MAG.* 3 ,E8.2,2A1) 
LARGEST AUX- 3 , II, 3 CORR. MAG.* \E8.2,2A1) 

AVE. AUX- 3 ,11, 3 CORRELATIOI* 3 ,2(1X,E11 .4) ,2A1) 
AVE. AUX- 3 ,11, 3 DC OFFSET* 3 ,2(1X,E10. 3) ,2A1) 

# SCABS FOR OFFSET* 10; t SCABS FOR AVE. CORR. 5 


M3 
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X ,2A1) 

911 FORMAT (' ARRAY OUT DC OFFSET= » ,2(1X,E10.3) P 2A1) 

950 TYPE*, 'I DISABLED THIS OPTIOH’ 

GOTO 25 

970 CALL RDPRAT 

IF (IFL .EQ. 1) GOTO 9999 
GOTO 25 
9999 IFL=0 

RETURI 
EID 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C Subroutine rdprat , to calculate and output the S/J,S/I ratios, 

C and simulated aidelobe levels associated vith a given exper- 
C linen t . 

SUBROUTIHE RDPRAT 

COHHOI/SPVRS/PVR , PISE , ITER, IFL , OFLAG 
IITEGER IFL, OFLAG 

REAL PWR(3,3) , PISE (3) ,SLL(3) , PRAT (3, 2) ,SIR(3) ,SIR 
PRAT(1 , 1 ) =10 ,0*ALOG10(PVR(1 , 1 ) /PVR(2 , 1 ) ) 
PRAT(l,2)*10.0*ALOG10(PVR(l,l)/PVR(3,l)) 

PRAT (2 ,1 )*10 .0*ALQG10(PVR(2 , 2) /PWR(1 , 2) ) 

PRAT (2 , 2)=*10 .0*AL0G10(PVR(2 , 2) /PWR(3, 2) ) 

PRAT (3 , 1 )=10.0*ALOG10(PWR(3,3) /PWR(1 ,3) ) 

PRAT (3 , 2) =10.0*ALOG10(PWR(3 ,3) /PWR( 2 ,3) ) 

SIR=10 .0*ALOG10(PWR(3,3)/(PWR(1 , 3)+PVR(2,3) ) ) 

SLL (1 ) ^AMIHl (PRATCl , 1 ) ,PRAT(1 ,2) ) 

SLL(2)*AMII1 (PRAT (2 ,1) , PRAT (2, 2) ) 

SLL(3)*AMII1 (PRAT (3,1 ) , PRAT (3, 2) ) 

C The factor -7.27 is due to the average poser of our pulsed 
C sinusoids, Pave=(Ac**2)/2* (duty cycle). The quantity PVR(I,Z) 

C represents (Ac**2)/2. The duty cycle cancels when taking signal to 
C signal ratios, but must be included for SIR computations. 

C UPDATE 8-9-88 BY DILSAVOR: THIS FACTOR IS TAKE! IITO ACCOUIT 

C BY SCALIIG PISE BY 3.0 AS SOOI AS IT IS CALCULATED II JVISC 
DO 18 1*1,3 

18 SHR(I)=10. 0*AL0G10(PWR(I ,I)/PISE(I) ) 

C18 SIR(I)=*10. 0*ALOG10(PVR(I ,I)/PISE(I) ) -7. 27 

C IF (OFLAG .IE. 1) GOTO 990 

TYPE 899 

IF (IFL .EQ. 1) TYPE 898 
49 DO 100 J=i , 2 

K=(-1)**J 
TYPE 900, J 

TYPE 901, J, J-K,PRAT(J,1) , J , PRAT( J, 2 ) , J , SIR(J) 

100 COITIIUE 

TYPE 902 

TYPE 903, PRATO, 1), PRATO, 2) ,S5R(3) 

TYPE 904, SIR 
RI1IRM=SIR(3) -PRATO, 1) 

RI2IRN=SIR(3)-PRAT(3,2) 

RLATT1=10.0*ALDG10( PWR(1 ,1)/PWR(1 ,3) ) 

RLATT2=10.0*AL0G10( PWR(2 , 2)/PWR(2 ,3) ) 

TYPE 905,RLATT1 ,RLATT2 
TYPE 908, RI1NRM ,RI2IRM 
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898 


GOTO 990 
FORMAT ( 1 

Bafor* Adaptation: 1 

,2A1) 

899 

900 


FORMAT (//) 

FORMAT (' AUXILIARY ELEMEIT » 

,11, > PARAHETERS; 

901 


FORMAT (' 


,F8.4,' dB ’ , ’ J',11 


X 

» dB V J 

r», 11, '/EOISE- ' ,F8.4 

, > dB\2Al) 

902 


FORM AT (' 

ARRAY OUTPUT (MAII AITEIIA IF WEI0HTS*0) 

903 

X 

,2A1) 
FORMAT ( 1 

D/Jl* \F8.4, 1 dB 

D/J2- '.F8.4,’ 

904 

X 

F8.4, » dB 1 
FORMAT (* 

',2*1) 

D/I-D/CJ1+J2)- 

F8.4 , ’ dB’,2Al) 

906 


FORMAT (' 

LEAKAGE ATTEIS: AUX1 

- H* ’ , Ffl . 2 , ’ AUX2-H«' 

906 


FORMAT (» 

I1IR IE HAIE-\F8.2, 

> I2IR IE HAIE-'.FB. 

990 


RETURI 

EID 




D/I 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c SUBPROGRAM TO APPLY I,Q WEIGHTS USIIG DAC-11 D/A COIVKRTERS 
C PROGRAHMIIG BOTH MAGIITUDE AID PHASE 

SUBROUTIIE JWWT 

COHMOI/VTBLR/LDAO.VI ,VQ, VMX ,VHY,IV 
REAL VI (2) ,VQ(2) , VMX (40 ,4) ,VMY (40,4) 

REAL A (4) , XLSB ,MAG, PHI 

IITEGER ID(4) , IV<4) , J , IADDRI, IADDRQ ,LDA0 

XLSB“20. /4095 
DO 150 J=l,2 
K=2*J 
L-R-l 

BO TI=VCV( VMX(1 ,L) ,VMY(1 ,L) ,IW(L) P ABS(WI(J))) 

TQ*WCV ( VMX ( 1 , K) ,VHY(1.K) ,IV(R) , ABS(VQ(J) )) 

VI**SIGI (TI , VI (J) ) 

VQ*SIGI(TQ ,WQ(J) ) 

C TYPE *,VI(J),VI,VQ(J),VQ 

80 IDI*IIT(VI/XLSB+"4000+0 . 5) 

IF(J .EQ. 1) LDAO*IDI 
IDQ*IIT(VQ/XLSB+”4000+0 . B) 

IADDRI*” 170440+ (J~l ) *4 
IADDRQ* I ADDRI+2 
CALL IPORE( IADDRI ,IDI) 

CALL IPOKE (IADDRQ, IDQ) 

150 COITIIUE 

099 RETURI 
EVD 

FUICTIOI WCV (X f Y ,I,Y1 ) 

REAL 1(40) ,Y(40) , Y1 
DO 100 J*1,I 
q-Y(J) 

R*Y (J+l) 

IF (Q.LE.Y1 .AID. Yl.LE.R) GOTO 1B0 

100 COITIIUE 

160 SLPE=(R-Q)/(X(J+1)-X(J)) 

VCV*(Y1-Q)/SLPE+X(J) 

RETURI 

EID 
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